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MODELLING ORGANIC COMPOUND REACTIVITY IN CYTOCHROME P450 MEDIATED REACTIONS 



RELATED APPLICATIONS 

This application claims benefit of priority from U.S. Provisional Application 
60/095,460, filed August 5, 1998, which is hereby incorporated by reference as if fully 
10 set forth. 



FIELD OF THE INVENTIQN 

The present invention is directed to the use of computational and other 

1 5 experimental data in the design of new pharmaceuticals, and in predicting the 

metabolism profiles thereof Computational and other information is used to further 
understand drug metabolism and toxicology, particularly in relation to 
monooxygenase enzymes, such as those of the cytochrome P450 system, that are 
involved in drug metabolism. Information derived according to the practice of the 

20 invention is useful in predicting the clearance or half-life of drugs, the propensity for 
drug interactions, and the nature and toxicity of byproducts resulting from their 
metabolism. The invention provides novel and powerful new approaches to drug 
design. 

25 BACKGROUND OF THE INVENTION 

The opportunity to improve the efficiency of the discoveiy phase of drug 
discovery is widely recognized in the pharmaceutical industry. Combinatorial 
chemistry and parallel synthesis methodologies, genomics, robotics, miniaturization, 
high-throughput screening and information technology together have stimulated an 

30 explosion of potential new lead compounds. However, limited expertise and 
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resources are available to move a compound from "candidate" to "lead". This is 

partially because the development of new technologies to address the later stages of 

the lead optimization process have not kept pace with the combinatorial technologies 

developed for synthesis and screening. As a result, there is a bottleneck in the drug 

5 discovery process that begins with lead optimization and extends all the way to the 

selection of clinical development candidates. 

A primary consideration in the area of lead development is a compoimd's 
metabolic fate. It would be of great value to the pharmaceutical industry if the 
discovery of such lead substances could be accelerated by design approaches that 

10 minimize reliance on screening synthesized compounds, but instead take advantage of 
quantifiable chemical or biochemical properties of a molecule in order to predict its 
metabolic characteristics. Predictive models can not only help alleviate the bottleneck 
in lead optimization but to facilitate the design and selection of drug candidates with 
not just adequate but with optimal absorption, distribution, metabolism and 

1 5 excretion/phannacokinetic (ADME/PK) profiles for progression to the later, more 
costly stages of the drug discovery/development process. Compounds thus selected 
will not only have a greater chance of development success but will ultimately lead to 
better medicines. 

In connection with the design of optimized pharmaceutical compounds, one 
20 particular area of interest involves a class of enzymes known as heme-containing 
monooxygenases, mixed fimction oxygenases, or alternatively cytochrome P-450s. 
These enzymes, which are abundant in the liver, act on a very broad range of substrate 
compounds, which is a trait very unusual for enzymes. These enzymes react with 
molecular oxygen such that one of the oxygen atoms is reduced to water, and the 
25 other is inserted into a substrate organic compound. This metabolic reaction may also 
be referred to as a hydroxylation reaction. 

The cytochrome P450 (CYP) enzymes comprise a superfamily of heme- 
containing enzymes that consists of more than 700 individual isoforms, and are found 

30 in plant, bacterial and mammalian species (Nelson et al. (1996) Pharmacogenetics 6, 
1-42). These enzymes fimction mainly as monooxygenases (Wislocki et al. (1980) in 
Enzymatic Basis of Detoxification (Jakoby, W. B., Ed.) pp 135-82, Academic, New 
York). In mammals, they are responsible for the metabolism of certain endogenous as 
well as exogenous compounds (Gonzalez, F. J. (1992) Trends Pharmacol Sci. 13 y 

35 346-52). The catalytic cycle for aliphatic hydroxylation by mammalian CYP is 
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described here briefly and shown in Figure 1 . Substrate binding changes the 
equilibrium of the heme iron from the low spin to the high spin state (step 1). This 
change lowers the reduction potential for the iron and facilitates electron transfer from 
NADPH via another enzyme, cytochrome P450 reductase (step 2). Molecular oxygen 
5 binds and is reduced by one electron as iron changes from the ferrous to the ferric 
state (steps 3 and 4). A second electron reduction of oxygen occurs and a peroxy 
intermediate is formed (step 5). The peroxy species undergoes heterolytic cleavage: 
one atom of oxygen leaves as a hydroxyl anion and the other forms a reactive oxygen 
species which is coordinated with the iron (step 6). Oxygen is transferred to the 

1 0 substrate (steps 7 and 8) and product dissociates from the enzyme (step 9). The first 
three steps of the cycle have been characterized spectroscopically. The next three 
steps occur rapidly, and have proven difficult to measure. At least two mechanisms 
have been proposed for the step(s) of oxygen transfer (depicted as steps 7 and 8 in 
Figure 1), and are described briefly below. 

1 5 The consensus mechanism for oxygen transfer is a nonconcerted reaction. In 

this context, nonconcerted implies that there are two distinct steps and that each step 
has its own transition state. There is evidence that step 7 (Figure 1) is abstraction of a 
hydrogen atom from the substrate by the reactive oxygen, which yields a carbon- 
based radical and an iron bound hydroxyl radical (White et al. (1980) Ann. Rev. of 

20 Biochem. 49, 3 15-56). The next step (step 8 in Figure 1) is rapid recombination of the 
two radical species, the "oxygen rebound" step. The magnitudes of isotope effects 
(Groves et al. (1978) Biochemical & Biophysical Research Communications 81, 154- 
60; and Hjelmeland et al. (1977) Biochem, Biophys, Res. Commun. 76, 541-9), loss of 
stereoselectivity (Groves et al.. Ibid; and White et al. (1986) J. Am. Chem. Soc. 108, 

25 6024-3 1), and evidence for rearrangement of the radical-like product of the first step 
(Groves et al. (1984) J. Am. Chem. Soc. 106, 2177-81) in various CYP-mediated 
reactions support the case for the oxygen rebound mechanism. 

Suitable substrates for CYP include steroids, prostaglandins, fatty acids, and 
exogenous drugs, pesticides and other toxic environmental contaminants including 

30 many carcinogens. Hydroxylation reactions are often the first step in the metabolism 
of foreign substances leading, for example, to the inactivation of administered 
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pharmaceuticals. Depending upon the mechanism of action of a drug, and its 
toxicological profile, it may be desirable to accelerate or delay its breakdown once it 
enters the body. Additionally other potential pharmaceuticals may be too toxic to 
administer, but appropriate structural modifications thereof may lead, for example, to 
5 structures that are decomposed to different, and less toxic, metabolites. 

For example, metabolism of phenylacetonitrile at the benzylic position (Figure 
2, see arrow) causes the release of the toxic metabolite cyanide, whereas aromatic 
oxidation leads to a less toxic product (Silver et al. (1982) Drug Metab. Dispos. 10, 
495-8). The metabolism of benzo(12)pyrene (Franchetti et al. (1995) J. Med, Chem. 

10 38, 3829-37) by CYP and epoxide hydrolase yields several metabolites including the 
extremely carcinogenic compound 7(R),8(S>dihydrodiol 9(S),10(R)-epoxide. These 
are examples where regioselectivity and stereoselectivity are important determinants 
of toxicity. In addition, potential drug interactions can be caused by 

differences in metabolism of multiple drugs by a single CYP isoform, as well as the 

1 5 induction or inhibition of individual CYP isoforms by drugs. Polymorphic enzyme 
expression among individuals may cause unusual patterns of drug metabolism, which 
can lead to imwanted side effects. 

As aforementioned, there is a tremendous need to enhance presently available 
methods for drug design to minimize dependence on the testing of randomly modified 

20 structures. But while progress has been made in the area of basic research, the field of 
predictive metabolism is not widespread in industry at this time. The two main 
determinants of enzymatic reactions are the steric and electronic characteristics of the 
enzyme and the substrate. Various electronic models have been developed that use 
computational (Grogan et al. (1992) Chem. Res. Toxicol 5, 548-52; Korzekwa et al. 

25 (1990) J. Am, Chem, Soc. 112, 7042-6; and Yin et al. (1995) Proa Natl. Acad, Set 
USA 92, 1 1076-80) and chemical (Karki et al. (1995) J. Am. Chem. Soc. 117, 3657- 
64; and Manchester et al. (1997) J. Am, Chem. Soc. 119, 5069-70) approaches for die 
successfiil prediction of relative reaction rates and isotope effect profiles for some 
classes of small compounds which are metabolized by CYP. The smsdl sizes of the 

30 compounds reduce, or nearly eliminate, the steric contribution to the outcome of 
metabolism. 
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In some cases steric factors are the main determinants of CYP— substrate 
interactions (Jones et al. (1995) Biochemistry 34, 6956-61). Computational methods, 
which include molecular dynamics and molecular modeling, are currently employed 
to probe the interactions between the substrate and the enzyme (Jones et al. (1993) J. 
5 Am, Chem, Soc. 115, 381-7). These methods rely on the crystal structure of a CYP 
protein, but only a few crystal structures exist for four (soluble) bacterial CYP 
isoforms (Ravichandran et al. (1993) Science 261, 731-6; Poulos et al. (1987) 
Journal of Molecular Biology 195, 687-700; Hasemann et al. (1994) Journal of 
Molecular Biology 236, 1 169-85; and Cupp-Vickery et al. (1995) Structural Biology 

10 2, 144-53). Crystal structures for the (membrane bound) mammalian isoforms are not 
available. Although there is low sequence homology between bacterial and 
mammalian CYPs, current data suggests that some degree of tertiary structure is 
conserved (Korzekwa et al. (1993) Pharmacogenetics 3, 1-18). With the bacterial 
isoforms as a template, homology modeling techniques have been used to generate 

1 5 structural models of the mammalian isoforms (Hasemann et al. (1 995) Structure J, 
41-62). The structural models have proven successful in locating residues near the 
active sites of mammalian CYPs (Szklarz et al. (1995) Biochemistry 34, 14312-22). 
The results have initiated site-directed mutagenesis studies, which have provided 
additional information about the structural details of the active site (Kobayashi et al. 

20 (1998) Biochemistry 37, 6679-88). This additional information has in turn been used 
to refine the homology models. The challenge remains to find a model that combines 
structural and electronic information for the prediction of the outcome of mammalian 
CYP-mediated reactions. 

The present invention is directed to the use of computational models, in 

25 concert with other experimental data, to identify probable therapeutic compounds 

having improved therapeutic or toxicological profiles. CYP enzymes are amenable to 
the development of electronically based or quantum chemical predictive 
computational models for two reasons. 1) Many isoforms display low substrate 
specificity due to weak enzyme-substrate interactions. 2) The catalytic step of the 

30 enzyme family is thought to be the heterolytic cleavage of molecular oxygen (step 6 
in Figure 1). Therefore, for substrates that have weak interactions with the binding 
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site of the enzyme, it is possible that the electronic features of the active oxygen and 
the inherent reactivities of various functional groups on a single molecule will dictate 
relative rates of metabolism. 

The invention thus provides a system whereby computational and 

5 experimental data can be combined to define a unified model of drug reactivity. The 
methods thus provide for computational models for predicting rates and 
regiospecificity of drug metabolism. Although the approach is described in reference 
to properties resultant from interaction with monooxygenase enzymes such as CYP, it 
will be appreciated that the general principles are applicable in a wide variety of 

1 0 contexts. 



SUMMARY OF THE INVENTION 

The present invention is directed to methods that use experimental data, and 
15 also quantum mechanical (includiag electronic configuration) data, to parameterize 
the predicted reactivities of various substances, and their metabolites and precursors, 
upon interaction with an enzyme. The methods are most applicable for enzymes with 
broad substrate specificity (or low substrate selectivity). Examples of such enzymes 
include monooxygenases, glucoronyl transferases, and glutathione transferases. 
20 Sample substances for which parameterization is desirable include various 

drug compounds or pharmaceutically active agents as well as any molecule 
introduced (such as by ingestion or inhalation) into a living orgaiusm. Upon such 
introduction, the substances may xmdergo reactions with various enzymes, including 
the monooxygenase cytochrome P-450 (CYP). Substrate reactions with CYP type 
25 enzymes, include for example, resultant hydrogen atom abstraction, aromatic 
oxidation, and metabolism at carbonyi groups or at heteroatoms. 

Thus one aspect of the invention involves analysis of regioselectivities of 
various possible enzymatic reactions, and of the imderlying chemical reactions, to 
identify structural features in pharmaceuticals that provide for advantages in non- 
30 adverse metabolism or bioactivation to toxic metabolites, half-life and dosing. As 
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surprisingly described below, the regioselectivity of the underlying chemical reactions 
is often of greatest importance in such determinations. 

In an uncatalyzed reaction, the regioselectivity of the reaction reflects the 
energetic differences for reaction of the various positions of the involved molecules, 
5 In an enzyme catalyzed reaction, however, the substrate compound is often bound to 
the enzyme in a very specific fashion, thereby changing the regioselectivity. The 
separation of reactivity differences intrinsic to the substrate from those attributed to 
steric interactions imposed by the enzyme is generally very difficult Stated 
differently, it would generally be difficult to predict the effect of a structural 

10 modification of a target compound upon its metabolism by ena^ymes, given that the 
molecular mechanisms of its metabolism involve an intimate combination of effects 
contributed by the metabolizing enzymes and the compounds themselves. 

According to the practice of the present invention however, it has been 
determined that, v^dth respect to metabolism of various compoimds (including 

1 5 potential pharmaceutical substances), large amoimts of predictive information 

pertaining to optimized metabolism by enzymes with broad substrate specificity are, 
in fact, directly derivable fi-om the regioselectivity of imderlying electronic structure 
of the reactants. As such, the information may be derived without requiring reference 
to the enzymes. This greatly facilitates computational analysis of optimized structural 

20 features in newly designed therapeutics. 

Without being limited as to theory, it is believed that certain enzymes, such as 
CYP, have evolved to metabolize a wide range of substances not otherwise readily 
metabolized or detoxified. As a result, and in order to provide a broad range of 
response, the normal high selectivity of enzymes for particular substances, and the 

25 steric requirements normally found at enzyme active sites, have been lessened, thus 
substantially enhancing the value of computational information that may be derived 
solely by reference to the electronic configuration of the metabolized substrates. 

Although the above concepts are primarily described herein with respect to 
human pharmaceuticals and CYP, the methods of the invention are also applicable to 

30 any enzyme with broad substrate specificity for the design of industrial or agricultural 
chemicals, and the like, where interest is focused not on optimizing dose or 
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phannacodynamics, but rather avoiding metabolism {bioactivation or inactivation) of 
substances that may come in contact with the human body (such as pesticides or 
pharmaceutical agents) where it is intended to avoid undesirable metabolism. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 depicts a proposed catalytic cycle for cytochrome P450. 

Figure 2 shows cytochrome P450 substrates which can be metabolized to toxic 
10 compounds. 

Figure 3 depicts results of hydrogen abstraction reactions with the methoxy radical. 
Figure 4 shows correction of AMI heats of reaction to match DFT values 

15 

Figure 5 shows experimental AAG versus predicted AAH for ratios of hydrogen 
abstraction and aromatic oxidation reactions. 

Figure 6 shows isotope effect profiles (lEPs) for six /?ar^ar-substituted toluenes with 
20 expressed CYPs. 

Figure 7 shows CYP substrates used for the measurement of regioselectivity. 

Figure 8 shows the formation of /?-cresol and /?-methoxyben2yl alcohol from p- 
25 methylanisole by P450. 

Figure 9 shows two metabolites from the cytochrome P450 mediated metabolism of 
j9-methylanisole. 

30 Figure 10 shows two metabolites from the cytochrome P450 mediated metabolism of 
a-chloromethyl-/?-xylene. 
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Figure 1 1 shows five cytochrome P450 substrates of varying size. 

Figure 12 shows resonance stabilization towards a carbon based radical for o- 
5 methylanisole. 

Figure 13 shows substrates used to parameterize aromatic oxidation. 

Figure 14 shows a correlation between experimental reaction rates with t-butoxy 
1 0 radical and DFT energies. 

DETAILED DESCRIPTION OF THE INVENTION 

As described above, the present invention is directed to methods of 

1 5 determining the electronic tendency for metabolism of a compoimd by enzymes with 
broad substrate specificity. Such methods use experimental data, and also quantum 
mechanical (including electronic configuration) data, to parameterize the predicted 
reactivities of various substances upon reaction with such enzymes. While such 
reactions may include steric factors in addition to electronic factors, or even be 

20 sterically driven, the methods of the invention provide a means of determining 
enzymatic metabolism of various substances. 

In one embodiment of the invention, such methods comprise the steps of 

a) determining, without reference to said enzyme, the regiospecific reactivity of 
each functional group in said compound towards metabolic transformation by any 

25 reaction(s) mediated by said enzyme; 

b) deriving one or more equations relating the activation energy for each enzyme 
mediated reaction to quantum chemical descriptors; 

c) determining, j&om said one or more equations, the metabolism of said 
compound by predicting the activation energies for each reaction; 

30 wherein said determined electronic tendency provides the electronic component 

for said reactions. 
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As such, the methods determine the metabolism of various compounds by use 
of the regiospecific reactivity of functional groups, such as C-H, C-C, C=C, CsC, 
C=0, C-N, C=N, -S-, -N-, -N= -CHO, -OH, and -COH, in combination vdth 
equations relating the activation energy for each enzyme mediated reaction to 
5 quantum chemical descriptors. The activation energy for each enzyme mediated 
reaction may further be related to the heat of reaction with said compound. One 
means for deriving the equations is by use of semiempirical methods. 

The equations can be refined by parameterization with experimental results to 
correct errors in activation energies due to the use of semiempirical methods. Such 
10 experimental parameterization may be accomplished by identifying the relative rates 
of substrate rotation in the active site of the enzyme by determining the isotope effect 
profiles for said reaction(s). 

The equations can also be refined by calculating high level quantum chemical 
parameters to correct the errors in activation energies. The use of high level quantum 
15 chemical parameters may be used in combination with the experimental 

parameterization described above, and includes the use of ab initio methods to correct 
the quantum chemical activation energy (AMI) and/or the spin contamination in said 
one or more equations. 

With such methods, it can be determined whether metabolism of a compound 
20 results in one or more of the following: 

(1) the production of an excess concentration of a toxic metabolite; 

(2) an excessive rate of metabolism of a pharmaceutically useful substance; 

(3) an inadequate rate of metabolism of a pharmaceutically useful substance; and 

(4) an excessive rate of metabolism of a non-pharmaceutically useful substance 
25 into a toxic metabolite. 

While the methods of the invention may be applied to any enzyme with broad 
substrate specificity, the followmg discussion is based on the use of a 
monooxygenase, CYP, as a representative example. As such, the predominant 
enzyme mediated reactions discussed concerns hydrogen atom abstraction and oxygen 
30 addition. These reactions include hydroxylation and aromatic oxidation. Obviously, 
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the use of other low substrate selectivity enzymes such as glucoronyl and glutathione 
transferases would focus on other chemical reactions. 

Essentially any new hydrophobic drug or xenobiotic will serve as a substrate 
for the cytochrome P450 (CYP) system. As such a predictive model for reactions 
5 catalyzed by these enzymes, in particular human CYP enzymes, would be an 

extremely useful tool. If the rate of a given CYP mediated oxidation of a substrate 
can be predicted, a number of properties could be anticipated, including: the 
clearance of drugs, the toxicity of xenobiotics, and the amounts of different 
metabolites that would come from a given substrate. These predictions will be useful 

10 in drug design, drug metabolism, and toxicology. 

Certain CYP enzyme characteristics render them suitable for theoretical 
calculations to be used to model their reactivity. First, the enzyme's principle 
catalytic steps involve the generation of an active oxygenating species. Other than by 
possibly initiating the steps, substrates are not directly involved. Upon generation of 

1 5 the active oxygenating species, the substrate is reacted in a manner that appears to be 
"chemical-like", such that only the broad regioselectivity of the reaction is influenced 
by the enzyme. Secondly, many CYP mediated reactions are unlikely to involve 
significant charge generation during the initial step of oxidation because the active 
oxygen is believed to have radical characteristics. Additionally, the CYP active site is 

20 thought to be a more hydrophobic than polar environment, so that gas-phase 
calculations would be a better approximation. 

With respect to the present invention, the following references are specifically 
incorporated by reference herein, as if fully set forth, and form a part of the present 
disclosure: J. Grogan et al., Chem. Res. Toxicol, Vol. 5, No. 4, pp. 548-552, 1992; H. 

25 Yin et al., Proc. Natl Acad Set USA, Vol. 92, No.5, pp. 1 1076-1 1080, 1995; J. Jones 
et al., "Predicting the Rates and Regiospecificity of Reactions Mediated by the P450 
Superfamily", inMs^W^* in Enzymology, Vol. 272, pp. 326-335, 1996; and K. 
Korzekwa et al., J. Am. Chem. Sac. , Vol 1 12, pp. 7042-7046, 1990. 

The predictive models for drug metabolism of the invention are based upon 

30 the application of quantum chemistry methods. These include semiempirical 
(electronic) methods and gaussian-based ah initio methods. These methods are 
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known in the art (see Lickowitz et al. ed., Reviews in Computational Chemistry II, 
VCH Publishers, 1991, pp. 313-315). To a first approximation, the electronic 
structure of a molecule is determined by the orientation of the atoms in space. Stable 
geometries of molecules are associated with atomic configurations that provide the 

5 lowest energy electronic configurations. Reactions occur through changes in atomic 
configurations from reactants to products. Along the reaction coordinate, the energy 
increases from the reactant geometry (stable, low energy configuration) to the 
transition state geometry (unstable, high energy configuration) and then decreases 
again to form products. 

1 0 Quantum chemical techniques attempt to model the electronic configurations 

and energies associated with atomic orientations. Approximate geometries can be 
optimized to stable geometries by minimization of the energy with respect to the 
atomic coordinates. Reactions can be modeled by transforming the reactant geometry 
to product geometry and minimizing all but one degree of freedom. 

1 5 Thus, the "heart" of a quantum chemical method is the method used to 

calculate the electronic structure for a given atomic configuration. The electronic 
configuration of a molecule is obtained by combining the atomic orbitals to form 
molecular orbitals. The equations for the electronic wavefimctions have been aroimd 
since the beginning of this century, but they are not amenable to solution. Therefore, 

20 different approximations are used in the solution of the equations. 

The complex part of the electronic structure calculation is associated with the 
multicenter electron repulsion terms. Semiempirical (electronic) methods 
approximate these terms with fimctions that are optimized to reproduce geometries 
and heats of formation. The resulting electronic configurations and energies are 

25 dependent on the initial parameterization of the semiempirical method. In contrast, ab 
initio methods make no use of experimental data. Instead these calculations 
approximate the atomic orbitals with combinations of gaussian fimctions that are 
amenable to integration. The more gaussian fimctions that are used (larger basis set), 
the more accurate the electronic representation. Also, the lower the energy, the better 

30 the calculation. 

Another important concept is that a lower energy can be obtained by mixing in 
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excited state electronic configurations. This is due to electron correlation, the 
tendency for electrons to keep away from each other. This energy gain is 
insignificant for most molecules, but are important for radicals, excited states and 
unstable geometries such as transition states. Accurate ab initio calculations require 
5 additional efforts to include correlation. 

Semiempirical calculations are fast, allowing very large molecules to be 
modeled. Ab initio calculations with large basis sets and correlation corrections are 
very slow, due primarily to the correlation correction. Calculations are presently 
limited to about ten heavy (non-hydrogen) atoms. 

10 An additional ab initio method based on density functional theory is also used 

in the invention. These methods include a correlation functional in the waveftmction 
calculation. Although these calculations are much faster than correlation-corrected- 
calculations, the exact form of the correlation functional is not implicit. Therefore, 
these functionals are parameterized, resulting in some of the same deficiencies of 

1 5 semiempirical calculations, i.e. a lower energy does not necessarily imply a better 

calculation. These calculations appear to be excellent for groimd states (i.e. reactants 
and products) and may be applicable for transition state calculations upon 
modification. 

The invention uses semiempirical calculations (AMI) to construct predictive 
20 models of drug metabolism. These models are based on the electronic tendency for 
metabolism, as determined by the reactions of substrates with small organic oxygen 
radicals. The speed of semiempirical calculations allows essentially any drug 
molecule to be modeled. Correlated (MP2) ab initio methods and density functional 
methods are used to provide more meaningful energies that can be used to validate 
25 and correct the semiempirical results. 

Based on the above, the present invention is directed to models for CYP 
mediated hydrogen atom abstraction and aromatic oxidation reactions based on 
reactions between a number of small radicals and substrates. The radicals used 
include p-nitrosophenoxy radical (PNR), t-butoxy radical, methoxy radical, HSO 
30 radical, and FO radical. Other radicals may be used for generating additional models 
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for hydrogen atom abstraction and aromatic oxidation reactions as well as other P450 
mediated reactions. 

PNR model limitations 
5 Our previously published model uses the p-nitrosophenoxy radical in 

hydrogen abstreiction reactions as a model for CYP mediated hydrogen atom 
abstraction. Using this model, the relative rates of reaction for several classes of 
small hydrophobic compounds can be accurately predicted. For example, the model 
has been able to predict the toxicity (LD50) of nitriles, the rates of in vitro metabolism 

1 0 of substituted toluenes and haloalkanes, and the rates of in vivo metabolism of 

halogenated hydrocarbon anesthetics. This model is based solely on computational 
data. Since only one kind of reaction is involved, parameterization of the linear 
model to generate absolute rates is not necessary, because the relationship betv^^een 
the relative rates will be maintained irrespective of the absolute values of the 

15 parameters. 

Although the PNR model is useful in predicting one class of CYP reactions 
(hydrogen atom abstractions), a complete predictive model for the CYP enzymes must 
also include other reaction types. These include aromatic oxidation, olefinic 
oxidation, carbonyl metabolism, N-oxidation and S-oxidation. Deficiencies with the 
20 previously published PNR model became apparent when other reactions were 
modeled. These problems are summarized as follows: 

1) Systematic errors within a computational method may not be consistent 
across different reaction types. Differences in systematic errors are expected since: 

a) The systematic errors are dependent on the nature of the potential 
25 energy surface. As an example, the absolute activation energies are greatly 

overestimated within the AMI formalism. The value of the computation is 
due to its ability to provide relative activation energies, 

b) Entropy terms are expected to cancel out within a series of related 
reactions, but will not cancel when different reactions are considered. 

30 2) Radical addition reactions have substantial spin contamination. This is a 

deficiency of the UHF (um-estricted Hartree-Fock) method, in which the wave 
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function for open shell systems are contaminated with higher spin state wave 
functions. This bestows more imcertainty on the relative energetics of the 
calculations. The transition states have <S^> values of approximately 2.1 and the 
tetrahedral intermediates have ,S^> values of approximately 1.1. Pure doublets have 

5 <S^ values of 0.75, indicating significant contamination with higher spin states. If 
the amount of spin contamination remains constant, as is the case for the substituted 
benzenes, the error is also likely to be constant and will not be problematic. However, 
calculations on polycyclic aromatic hydrocarbons show that the degree of spin 
contamination depends on the extent of conjugation of the radical, i.e., the degree of 

1 0 spin contamination is not size consistent. Therefore, different substrates will have 
different errors due to spin contamination. 

3) Although addition of the nitrosophenoxy radical to an aromatic carbon 
progresses smoothly through a transition state and on to a tetrahedral intermediate, the 
electronic state of the tetrahedral intermediate is not the ground state. Thus, the 

1 5 reactions from reactants to products and products to reactants occur on two different 
potential energy surfaces. In addition to complicating any transition state studies, the 
value of any Bronsted correlations will be in doubt, since the transition states are not 
necessarily intermediate between products and reactants. 

Therefore, it was necessary to use a completely different approach to develop 

20 more comprehensive predictive models for the CYP enzymes. This invention 

describes a method that uses experimental data and additional quantum chemical 
methods to parameterize models that will predict several types of CYP mediated 
oxidation reactions, includmg hydrogen atom abstraction, aromatic oxidation and the 
metabolism of carbonyl compounds. This model does not require the calculation of 

25 transition states for a model oxygen radical such as the nitrophenoxy radical. Instead, 
the model is parameterized using experimental regioselectivities of enzymatic 
reactions. Additional input for model development is provided by rates of reactions 
of chemical systems, and ab initio quantum chemical calculations using alkoxy 
radicals. 

30 

Parameterization with Enzymatic Systems 
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Based on a number of C YP characteristics described in the Examples below, 
the validity of deriving predictive information from the underlying electronic structure 
of the reactants is demonstrated. The characteristics include the demonstration by 
kinetic isotope effects (KJEs) that the active CYP oxygen is equivalent in a number of 
5 CYP isoforms; and the indication by intramolecular and intermolecular KIEs as well 
as the ratio of products from competing CYP mediated reactions on a model substrate 
that rapid substrate rotation permits use of electronic factors alone to predict 
regioselectivity. The latter point is of particular significance because it permits the 
use of model radicals in place of CYP enzymes in obtaining model data for deriving 

10 predictive models. 

Given these results, models for hydrogen atom abstraction and aromatic 
oxidation were derived upon application of quantum chemistry methods. 

Use of semiempirical (electronic) methods, however, result in incorrect 
absolute activation energy (AMI) barriers, making the methods inappropriate for 

1 5 comparing the energetics of different reaction mechanisms. But it is preferable to use 
models based on these methods to minimize the computational costs of the 
calculations. In order to calculate ratios of metabolites from a single compound and 
to include other reaction mechanisms, the AMI barriers need to be corrected by 
correlation with experimental data. Experimental data is obtained from metabolite 

20 ratios when a substrate xmdergoes rapid exchange between positions of metabolism. 
This rapid exchange, which is also referred to as rapid equilibrium, reflects the rapid 
substrate rotation in CYP mediated reactions and may be determined by KIE studies. 



Hydrogen atom abstraction 

25 Based on the above, an initial model for hydrogen atom abstraction by CYP of 

the invention was formulated based on the semiempirical quantum chemical 
activation energy (AMI) from quantum chemical descriptors of the substrate and the 
radical product. This can be expressed by the following formula of Equation 1, 
relating the activation energy (AHact) to the respective heat of reaction (AHr) for each 

30 hydrogen atom abstraction reaction and the ionization potential of the resulting radical 
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(IPrad). This allows for the prediction of AMI activation energies with a rqpid 
calculation of ground state properties. 

AHact = 2.60 + 0.22 (AHr) + 2.38 (IPrad) Eq- 1 

5 

This equation, when used alone, can reproduce the AMI activation energies, 
but the AMI formalism is known to overestimate these barriers. Thus, as described 
above, experimental data can be correlated with this equation to correct the 
overestimation. This is possible upon correlations with regression of experimental 

10 data on parameters that include functions of this AMI activation energy. The 

overestimated barriers are corrected by a proportionality constant of the regression 
equation. The improvement based on correlation with the regression is dependent on 
the number of experimental data points used in the regression. As such, the accuracy 
of the equation can be repeatedly improved by correlating with regressions based on 

15 increasing amounts of experimental data. 

Five substrates shown to undergo rapid equilibrium between positions of 
metabolism have been used to prepare data for regression and then correlation with 
Eq. 1 . The differences in activation energy calculated by the experimental 
regioselectivity of the substrates, along with the differences in AMI activation 

20 enthalpies, are given in Table 1 . Regression of the experimental data gives a slope 
correction for AMI activation energies of 0.75 ± 0.20. 



Table 1. Experimental and AMI energy differences. 



Substrate 


AAG (experimental) 


AAH(AM1,PNR) 


4-methylanisole 


1.1 


2.2 


2-methylanisole 


0.3 


2.3 


a-chloro-p-xylene 


-1 


-1.3 


Octane 


-3.13 


-2.6 


Valproate 


-2,48 


-2.6 
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Another source for correcting AMI activation energies is with high level ab 
initio calculations, which may be combined with experimental corrections. As stated 
above, t-butoxy and methoxy radicals can be used as models for CYP mediated 
hydrogen atom abstraction reactions. Table 2 gives the ab initio MP2/631 1+G* 
5 activation energies for methoxy radical mediated hydrogen atom abstraction, 

compared to the AMI results with the p-nitrosophenoxy radical. MP2 and 63 1 l+G* 
are ab initio quantum chemical formalisms known in the art. 



Table 2. Ab initio and AMI Activation Energies. 



Substrate 


AHmp2(MP2/6311+G*) 


AAHami 


Propene 


10.60318 


21.53 


methylamine 


9.078607 


17,78 


acetonitrile 


15.72245 


26.33 


ethane 


13.34009 


24.57 


methane 


16.46698 


26.99 



10 



Regression of this data gives Eq. 2, with an r^ = 0.95. 

AHmp2 = -6. 1 1 + 0.82 (AHami) Eq. 2 

15 The slope correction of 0.82 is similar to the experimental correction of 0.75 

described above. The intercept correction of -6.1 1 is essentially cosmetic because it 
cancels when relative free energy differences are calculated. Using the experimental 
correction of 0.75 and an intercept correction of -6.1 1 as corrections to Eq. 1 results in 
the experimentally corrected equation for hydrogen atom abstraction (Eq. 3). 

20 

AHact = -4. 1 6 + 0, 1 65 (AH)r +1.78 (IPrad) Eq. 3 

This equation gives the best fit of AMI heats of reaction and ionization 
potential to the experimental data in Table 1. It also brings the absolute values in line 
25 with the calculated MP2/63 1 1+G* activation energies. 



18 



BNSDOCID: <WO _ OOOSSOSAI I > 



wo 00/08205 



PCT/US99/177I3 



Given the above, predictions for the electronic tendency for hydrogen atom 
abstraction of a given compound of interest by C YP can be made by simply 
calculating the radical and substrate heats of formation, which provide the heat of 
reaction by subtracting the former from the latter, and the ionization potential of the 
5 radical. For substrates that rotate rapidly in the P450 active site, the model will 
predict regioselectivity. If steric effects are involved, the model will predict the 
electronic component of the reaction. 

The validity of this method to predict hydrogen atom abstraction is Observed 
for the CYP3 A4 enzyme. This enzyme has the broadest substrate specificity and 
1 0 apparentiy minimal steric interactions for most substrates. A representative list of 
CYP3A4 substrates is shown in Table 3. 



Table 3. Use of CYP3A4 Regioselectivity to Validate the Electronic Model 



15 



for Hydrogen Atom Abstraction. 



Compound 



Buprenorphine 

Cocaine 

Estrone 

Nifedipine 

Nitrendipine 

Morphine 

Ethylmorphine 

Dextronnethorphan 

Cyclobenzaprine 

Codeine 

Clozapine 

Tamoxifen 

Testosterone 

Salmetrol 

Quinidine 



Primary metabolism 

predicted to be 1" 

Electronic 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 



Primary metabolism Not in top 3 

predicted to be in top 3 

< > Steric 



X 
X 



In this Table, the model accurately predicted the primary metabolite for 12 of 
15 substrates. This suggesets that the oxidation of these substrates is electronically 
controlled. Compounds that are metabolized in deactivated positions are likely to be 
sterically controlled. 
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Aromatic Oxidation 

Aromatic oxidation is likely to occur through a radical addition mechanism. 
Unfortunately, addition of the p-nitrosophenoxy radical to an aromatic carbon 
progresses smoothly through a transition state and onto a tetrahedral intermediate 
5 where the electronic state of the intermediate is not the ground state. Thus, the 

reactions from reactants to products and products to reactants occur on two different 
potential energy surfaces. In addition to complicating transition state studies, the 
value of any Bronsted correlations for addition to substituted benzenes will be 
doubtful because the transition tates are not necessarily intermediate between products 

10 and reactants. 

Thus the p-nitrosophenoxy radical alone is not likely to serve as a versatile 
model for other CYP mediated oxidations. The difficulties in developing a quantum 
chemical model with the PNR model can be summarized as 1) inaccurate energies 
associated with AMI calculations, 2) spin contamination for unrestricted Hartree- 

1 5 Fock (UHF) calculations, and 3) surface crossings for the p-nitrosophenoxy radical 
model. These problems have been solved by 1) using experimental parameterization, 
and 2) use of high level ab initio calculations to guide the experimental 
parameterization, and 3) use of methoxy radical as the model oxidant.. 

The validity of using experimental parameterization is presented by the 

20 following, which shows the likely rapid rotation of substrates undergoing CYP 
mediated aromatic oxidation and the validity of using radicals to model CYP 
mediated aromatic oxidation reactions. 

To address the issue of inaccurate energies associated with AMI calculations, 
small oxygen radicals other than the p-nitrosophenoxy radical were tested within the 

25 AMI fomialism. HSO and FO radicals were foimd to add smoothly to aromatic 
compounds using the AMI fomialism. In contrast to the p-nitrosophenoxy radical 
previously used, HSO and FO radicals remain on a single potential energy surface. 
Both radicals provide excellent modified Bronsted correlations for addition to 
substituted benzenes, with values of 0.98 for HSO and 0.99 for FO. These 

30 excellent correlations for both the meta and para positions indicate that models can be 
developed to describe diverse electronic characteristics. This data also suggest that 
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the activation energies will simply be a function of heats of reaction and do not need 
to be corrected for resonance. 

The issue of significant amount of spin contamination along the reaction 
coordinate for UHF calculations is a more difficult problem. This reflects a 
5 deficiency of the UHF method, in which the wave fimction for open shell systems are 
contaminated with higher spin state wave functions. This results in added imcertainty 
in the relative energetics of the calculations. The transition states have <S^ values of 
approximately 1 .2 and the tetrahedral intermediates have <S^ values of 
approximately 1 . L Pure doublets have <S^ values of 0.75, indicating significant 

1 0 contamination with higher spin states. If the amount of spin contamination remains 
constant, as is the case for substituted benzenes, the error is also likely to be constant 
and not problematic. Calculations on polycyclic aromatic hydrocarbons, however, 
show that the degree of spin contamination depends on the extent of conjugation of 
the radical, i.e. the degree of spin contamination is not size consistent. Thus 

1 5 correction for spin contamination is necessary. 

One means of correcting for spin contamination is with ab initio calculations. 
MP2, MP4, and density functional ab initio calculations have been used to obtain 
ground state substrates, and when possible and necessary, transition states. Obtaining 
reasonably accurate energetics requires at least an MP2 level of theory with a 6-3 IG* 

20 basis set. MP4, 6-3 IG*, and DFT are ab initio quantum chemical formalisms known 
in the art. Using this basis set, hydrogen atoms were abstracted fi*om a series of 
substrates using /er/-butoxy and methoxy radicals. The methoxy radical adds 
smoothly to aromatic systems. The MP2 calculations have minimal spin 
contamination after projection. 

25 The modified Bronsted correlations for the methoxy radical calculation are 

shown in Figure 3. Similar results were obtained for the /er/-butoxy radical 
calculations. The methoxy radical was used since it is substantially smaller than the 
/er/-butoxy radical. Since the upper limit for the number of non-hydrogen atoms in a 
molecule is approximately ten, many more reactions can be modeled if the methoxy 

30 radical is used. The energetics for these two molecules are likely to be similar since 
the experimental OH bond dissociation energies are almost identical. As with the 
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semi-empirical PNR model, reasonable modified Bronsted correlations are obtained 
for both methoxy (R^=0.91) and tert-butoxy radicals (R^=0.92). Therefore, methoxy 
radical calculations (6-31G*/MP2) may provide a reasonable model for P450 
oxidations- These results suggest that ab initio calculations with alkoxy radicals can 
5 be used to guide the parameterization of the computational models. 

Another method of correcting for spin contamination is with density functional 
calculations. These methods are fast, high level ab initio methods that can provide 
correlated results with very litde spin contamination. Although their value in 
transition state calculations is uncertain, they can provide good uncontaminated 

10 ground state energies. For example, the bond dissociation energies for seven carbon 
hydrogen bonds were found proportional to the spin projected MP2 bond dissociation 
energies (R^=0.93), and an excellent correlation between /er/-butoxy radical 
abstraction rates and DFT calculations has been obtained (see below). DPT 
calculations have no spin contamination. 

15 Therefore, an approach is to use density functional ab initio methods to correct 

for spin contamination in heats of reaction. These corrected heats of reaction can then 
be used to predict activation energies for aromatic oxidation. Figure 4 shows the 
correlation between density functional (B3LYP) and AM 1 heats of reaction for 
methoxy radical addition to aromatic rings. The parameters for the fit is given in Eq. 

20 4, where a value AS^ is a correction for UHF spin contamination in the AMI 
calculation. 

AHr, DFT = 19.2 + 1.32 (AHr, ami) + 5.47 (AS^) Eq. 4 

25 The generation of an accurate combined model for hydrogen atom abstraction 

and aromatic oxidation requires consideration of competing oxidation reactions 
hydrogen atom abstraction reactions since. The relationship between these competing 
reactions can be expressed by Equations 5 and 6, where a is the intercept correction 
for aromatic oxidation, b is the slope term for aromatic oxidation, the term a(l-y) is 

30 used with y=l when an intercept correction is not needed (as in Eq. 5) and y=0 when 
an intercept correction is needed (as in Eq.6). These equations permit simultaneous 
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solution of data sets containing aromatic oxidation with and without hydrogen 
abstraction ratios. 

In 5^ = a(l-y) + b(AAHR.a«)m) Eq. 5 



In ?2I2E = a(l-y) + b(AHR,a„>m)-AHHabs Eq. 6 

PHabs 

1 0 These equations may be solved with the data in Table 4, which provides 

experimental data for aromatic oxidation reactions with a number of substrates. 



Table 4. Experimental Product Ratios for Aromatic Oxidation 



Substrate 


Ratio 
Description 


Experimental 
Ratio 


AAH(AM1) 

Habs 


AHr(AM1) 
aromatic 


AAHr (AMI) 
aromatic 


Toluene 


Arom/benz 


0.13 


12.32 


-7.19 




Anisole 


Arom/benz 


10.6 


11.29 


-6.39 




Ethylbenzene 


Arom/benz 


0.086 


10.21 


-6.83 




2-methylanisole 


Arom/benz 


1.76 


10.52 


-7.12 




Phenanthrene 


I-OHP/9-OHP 


0.61 






-1.90 


Phenanthrene 


4-OHP/9-OHP 


0.35 






-3.40 



1 5 Upon solution, the data in Table 4 gives parameters of a=ll .72 and b=0.28. 

Solving for aromatic oxidation on the corrected hydrogen abstraction scale results in 
Eq. 7. 

AHan^n. = 11 .25 + 0.24 (AHr. arom) Eq. 7 
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Combining Equations 7 and 4 results in Eq. 8, which predicts aromatic 
oxidation with AMI parameters. 

AHarom = 15.86 + 0.32 (AHr,ami) + 1.31 (AS^) Eq. 8 

5 

As with the equations for hydrogen abstraction provided above, the accuracy 
of the above equations can be repeatedly improved by use of increasing amounts of 
experimental data. 

10 Given the above, and similar to the situation for hydrogen atom abstraction, 

predictions for aromatic oxidation of a given compound of interest by C YP can be 
made by calculating the AMI heats of reaction and measuring the spin contamination 
of the resulting radical. Eq. 8 above, is used to predict the activation energy of the 
reaction. This activation energy can be directly compared with hydrogen abstraction 

15 reactions calculated by Eq. 3. Figure 5 shows the ability of the above methods to 
predict both hydrogen abstraction and aromatic oxidation. The energy differences 
obtained from both experimental and predicted data demonstrates the consistency of 
the predictive model. 

The equations correct models derived from pure AMI quantum chemical 

20 calculations to better match experimental data. Availability of additional 

experimental data will provide a means for further refinement of the coefficients in 
each equation. 

Other C YP mediated reactions 

25 CYP is capable of mediating additional reactions beyond hydrogen atom 

abstraction and aromatic oxidation. These reactions include olefinic oxidation, 
carbonyl metabolism, N-oxidation and S-oxidation, The formulation of predictive 
models for these reactions may be conducted in the same maimer as that provided 
above for hydrogen atom abstraction and aromatic oxidation. The mechanism of each 

30 reaction is first modeled and then separately parameterized by experimentation as 
provided above. 
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For example, carbonyl metabolism is likely to be hydrogen atom abstraction 
based. As such, rapidly rotating substrates capable of the carbonyl metabolism 
reactions of compounds of interest will be modeled with semiempirical (electronic) 
methods and then parameterized with ab initio methods. Similarly, N- and S- 
5 oxidation probably proceed by oxygen addition, which may also be modeled by 

rapidly rotating substrates undergoing oxygen addition. These substrates may then be 
used for experimental parameterization by both semiempirical and ab initio methods. 

Having described the invention with regard to C YP enzymes in general, it 
10 should be clear that all of the disclosed invention is applicable to human CYP 
enzymes, including, but not limited to CYP2E1, CYP3A4, CYP2B6, CYP2C8, 
CYP2C9, CYPl Al, CYP1A2, CYP2C19, CYP2D6, CYPIBI, and CYP2A6. 

The following examples are intended to support and illustrate but not to limit 

the invention. 

15 

Example 1 

General computing methods: 

A variety of software packages that are commercially available or otherwise in 
the public domain are useful in performing the quantum chemical calculations used in 

20 the practice of the invention. 

Both ab initio and semiempirical calculations can be performed using 
Gaussian 94 or Gaussian 98, which are available from Gaussian, Inc., Carnegie Office 
Park, Building 6, Suite 230, Carnegie, PA 15106, USA, Semiempirical and density 
functional calculations can be performed from within the Insight II software offered 

25 by Molecular Simulations Inc., 9685 Scranton Road, San Diego, CA, USA. 

Semiempirical calculations can also be performed using the conmaercial and public 
domain versions of MOP AC which is available from the Quantum Chemical Program 
Exchange, Indiana University, Creative Arts Building, Bloomington, Indiana, USA 
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Example 2 

Similarity in energetics of active oxygen in CYP enzymes: 

Due to heterogeneity in the CYP enzyme family, it was necessary to confirm 
5 that the energetics of any given reaction would be similar for each CYP. One means 
of confirming this is by comparing isotope effect profiles (lEPs) of a number of CYP 
enzymes to verify that the active oxygen in each energetically equivalent. For CYP 
mediated reactions, the physical characteristics of a compound and the nature of the 
active oxygen species dictate the magnitude of an intrinsic isotope effect. An lEP 

1 0 consists of isotope effect measurements for a series of structurally related compoimds 
that most nearly approximate intrinsic isotope effects using multiple CYPs. When 
intrinsic kinetic isotope effects (KIEs) are available, comparable trends in lEPs 
between different systems or different isoforms support similar mechanisms (Karki et 
al. (1995) J. Am. Chem. Soc. 117, 3657-64; Manchester et al. (1997) J. Am. Chem. 

15 Soc. 119, 5069-70; and Hemies et al. (1982) Biochemistry 21, 5106-1428). 

Isotope effects are reaction rate changes brought about by isotopic 
substitution. These usually carry mechanistic information. This invention is based on 
the comparison of relative reaction rates and enzyme kinetic constants for metabolism 
of CYP substrates that contain a carbon-hydrogen (C-H) bond to the metabolism of 

20 substrates that contain a carbon-deuterium (C-D) bond. A primary isotope effect 
refers to the rate change that is measured when the bond to the isotopic atom is 
broken- A secondary isotope effect refers to the rate change that is measured when 
the bond to the isotopic atom is not broken. 

Deuterium isotope effects arise from mass differences between hydrogen (H) 

25 and deuterium (D). The frequencies (v) for the stretching of C— H and C— D bonds are 
around 9 and 6.3 (xlO'^ s\ respectively (Isaacs, N. S. (1987) Physical Organic 
Chemistry, John Wiley & Sons, Inc., New York). The zero point energy of a bond, 
which is the vibrational energy remaining at 0 K, is defined by Ey = V2 hv, when V 
(the vibrational quantum number) = 0, and where h ~ Planck's constant. The zero 

30 point energy for a C-D bond is lower than that for a C— H bond. The isotopic 
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10 



15 



20 



difference in zero point energy decreases as V increases, and approaches zero in the 
transition state when bond dissociation is about to occur. Therefore, more energy is 
required to break a C-D bond than a C-H bond. 

In order to obtain mechanistic information from enzyme kinetic studies, 
intramolecular KIEs can be used. While this type of experimental design cannot be 
used for all enzymes, it is appropriate in cases where it can be used to provide a closer 
estimate of the intrinsic isotope effect, since the potential for masking is decreased. 
An intrinsic isotope effect is the maximimi isotope effect for a compound, and is 
dependent on only the inherent physical properties of the C-D and C-H bonds 
(Northrop, D. B. (1975) Biochemistry 14, 2644-51). Rate limiting steps in an enzyme 
catalytic cycle that occur before the isotopically sensitive step can have the effect of 
reducing, or masking, the true (or intrinsic) rate difference associated with isotopic 
substitution. 

For the intramolecular isotope effect experiments herein, a substrate is chosen 
that contains at least two symmetrical sites, with at least one site, but not all, 
isotopically labeled. A single concentration of substrate is incubated with enayme, 
and the ratio of labeled to unlabeled compound in the isotopic mixture of products (or 
the conresponding rates of formation of products) are used to calculate a KIE. 



Scheme I illustrates a simplified mechanism for an intramolecular isotope 

effect. 



Scheme 1 
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(Eq. 10) 



observed 




10 



15 



20 



Rate constants kj and k2 represent substrate binding to and debinding from 
enzyme (see step 1 in the Figure 1 catalytic cycle). Rate constant is shown as an 
irreversible step that represents steps that occur after substrate binding, up to and 
including oxygen activation (steps 2-6 in the catalytic cycle. Figure 1). An 'activated 
complex' between the en2yme and the substrate is formed (ES*) as a result of oxygen 
activation. Two possible orientations exist for the substrate in the active site 
(occurring between steps 6 and 7 of the catalytic cycle, Figiu-e 1). The orientations 
differ in the proximity of the H or D to the activated oxygen. Exchange between the 
positions is represented by the steps fcj, A4', and fe', and symbolizes rotation about 
the methyl group for a substrate of the form R-CHjD. Product release is represented 
by kg (step 9 in the catalytic cycle of Figure 1). 

The observed isotope effect is the experimentally measured ratio of Ph to Pd, 
and is described by Eq 9. The concentrations of the ES* species can be expressed in 
terms of rate constants to yield Eq 10. The intrinsic isotope effect, A;7h/^7d, and the 
ratio of kya to ^ determine the observed isotope effect. Equation 10 reveals that as 
the magnitude of A4 increases, the observed KIE will approach the intrinsic KIE. 
Since rotation about a methyl group is known to be fast, A4 is large relative to km» and 
the observed intramolecular KIE closely estimates the intrinsic KIE. 

The 'product' of the hydrogen atom abstraction step carried out by the iron- 
oxygen species [FeO]^^ in CYP can be thought of as a benzylic carbon radical, which 
allows one to evaluate the reaction coordinate for the hydrogen atom abstraction step. 
According to the Hammond postulate (Hammond, G. S. (1955) J. Am. Chem, Soc. 77^ 
334-40), a highly exothermic reaction will be expected to have a transition state 
resembling reactants, while an endothermic reaction will have a transition state more 
similar to products. As a substituent becomes more electron withdrawing, a benzylic 
radical is less stabilized. Consequently, the reaction is less exothermic (more 
endothermic), and the magnitude of the intrinsic isotope effect increases. 
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A three-center transition state complex and Melander and Westheimer's 
theory for the explanation of the magnitude of isotope effects (Melander, L. (1960) 
Isotope Effects on Reaction Rates, Ronald Press. New York; Westheimer, F. H. 
(1961) Chem. Rev, 61, 265-73; and Kresge, A. J. (1977) in Isotope Effects on Enzyme- 
5 Catalyzed Reactions (Cleland et al., Eds.) pp 37-63, University Park Press, Baltimore) 
can be used in the mechanistic interpretation of lEPs. The main concept of the theory 
is that the symmetric stretching vibration can supply isotope dependent zero-point 
energy to the transition-state, which can offset the isotope dependent zero-point 
energy of the ground state molecules. In a symmetrical transition-state, the force 

10 constants of the two partial bonds are equal, and the central atom is stationary. 

Therefore, replacing H with D will not change the moving mass, and vh will be the 
same as vd. The isotope dependent zero-point energy difference in the transition-state 
will be smallest for an isothermic (symmetric) reaction, leading to the largest isotope 
effect. When f\ ^fi, vh vd (i.e., for asymmetric transition states), and isotope effects 

1 5 will be smaller than the KIE for a symmetric reaction (Kresge, Ibid\ and Cleland, W. 
W. (1982) Methods Enzymol 87, 625-41). 

The CYP mediated primary KIE for octane is near the maximal possible value 
for hydrogen atom abstraction, therefore, a symmetrical transition state and an 
isothermic reaction are inferred. Thus, increasingly endothermic reactions relative to 

20 octane will yield increasingly smaller intrinsic ICIE, and increasingly exothermic 

reactions will yield increasingly smaller intrinsic KIE. A plot of KIE as a function of 
the extent of hydrogen transfer in the transition state would be a parabola. 

An lEP consists of isotope effect measurements in a system {e.g., chemical or 
enzymatic) using a series of compounds that are structurally related. To determine if 

25 the mechanisms of two reactions are similar, lEPs are preferred over a single 

measurement since the trend in isotope effects provides additional information about 
the energetics of the reaction and makes fortuitous agreement between the two 
systems unlikely (Karki et al. (1995) Xenobiotica 25, 71 1-24). When intrinsic KIEs 
are available, comparable trends in KIE profiles between different systems argue for 

30 similar mechanisms, while different trends are indicative of dissimilar mechanisms. 
As described above, maximum isotope effects will be observed with reactions that 
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have symmetrical reaction coordinates. Mechanisms with non-symmetrical reaction 
coordinates will yield smaller isotope effects. From this information it can be inferred 
that intrinsic KIEs that are measured for a series of compounds that differ in some 
systematic (e.g,, chemical or physical) way can serve as a sensitive probe of the 
5 energetics of the isotopically sensitive step. 

The lEPs for multiple CYP isoforms were measured in order to test the 
hypothesis that the active oxygen of multiple CYPs is conserved. The intramolecular 
isotope effects for six jt^am-substituted toluenes (p-X-Ph-CHiD: X = OCH3, CH2D, H, 
CI, Br, CN; Ph = phenyl) using expressed CYPs 1 A2, 2B1, 2E1 and 2C9 and purified 

10 P450cam are listed in Table 5. These experiments were performed with expressed 
C YP3 A4 also, but product formation was too low to accurately measure. Isotope 
effects measured with substrates that contain hydrogen and deuterium bonded to the 
same carbon atom permit close measures of the intrinsic isotope effects (see prior 
discussion of intramolecular isotope effects). Table 5 shows that for each isoform, the 

1 5 lEPs for the toluenes are practically identical: the values range from 4-11, and the 
magnitudes of the isotope effects have the same rank order. In addition, for each 
substrate the KIE for all five CYPs are in close agreement, which exemplifies the 
conserved nature of the active oxygen species. The data in Figure 6 illustrates the 
close agreement of intramolecular isotope effects for each substrate in multiple CYPs, 
20 which exemplifies the conserved nature of the active oxygen species. 



Table 5. Intramolecular Isotope Effects for Six Para-Substituted Toluenes 
Using Various CYP Isoforms 



/7-substituent 


2E1 


2B1 


1A2 


2C9 


P450cam 


OCH3 


4.24 (1) 


3,69 (3) 


4.64 (4) 


4.3 (1) 


4.44 (1) 


CHjD^ 


5.45 (5) 


6.23 (5) 


5.59 (1) 


5.9(1) 


6.0 (2) 


H 


6.1(1) 


7(1) 


6.1 (4) 


nd^ 


nd 


CI 


6.75 (2) 


8.1 (1) 


7.06 (2) 


6.2(1) 


6.5 (8) 


Br 


6.75 (3) 


8.02 (3) 


6.83 (8) 


6.9(2) 


8.3 (1) 


CN 


10.1 (1) 


11.9(3) 


10.1 (2) 


11.1(6) 


11.6(5) 



Each result is the average of three determinations and the numbers in 



30 



..0008?O5A1 J_> 



wo 00/08205 



PCTAJS99/17713 



parentheses indicate the standard deviation in the last significant digit of the mean 
isotope effect value. "/^-Xylene-a-^Hi-a'-^Hi was used as the substrate. * Not 
determined due to low levels of product formation. 

The lEPs herein provide evidence for a common reaction mechanism for 
aliphatic hydroxylation at the benzylic position in four mammalian CYP isoforms as 
well as one bacterial isoform. Inclusion of P450cam into the series of lEPs is the first 
extension of this application to a P450 outside of the Class II P450s. P450cam 
belongs to Class I, which consists of P450s that interact with two electron-transfer 
partners. Class II members differ in that they have only one protein as a reducing 
partner (Peterson et al. (1995) in Cytochrome P450: Structure, Mechanism, and 
Biochemistry (Ortiz de Montellano, P. R., Ed.) pp 151-80, Plenum Press, New York). 
The finding that the KIE profile for P450cam is not different fi-om the C VPs that 
belong to Class II is significant because it indicates that the mechanism of hydrogen 
atom abstraction is independent of the steps that occur prior to the irreversible 
formation of an active oxygen species. 

These results agree with other lEP data in the literature for different groups of 
compounds that are metabolized by CYP. Isotope effects for a series of N,N- 
dimethylanilines were measured using mammalian CYPs 2B1 (purified), 4B1 and 
1A2 (expressed), rat liver microsomes, and purified bacterial CYP102, and compared 
them to KIE in two different chemical systems (Karki et al. (1995) J. Am. Chem. Sac. 
117, 3657-6416). Practically identical lEPs for all of the P450s and the /er/-butoxyI 
radical system provide strong evidence for a similar mechanism of iV-dealkylation, 
i,e., initial hydrogen atom abstraction. Furthermore, the mechanism was conserved 
among all the CYPs that were tested, which includes both mammalian and bacterial 
members of the Class 11 P450s. Isotope effect profiles reported for various CYPs and 
for /er/-butoxy radical in several small molecules also support the mechanism of 
(initial) hydrogen atom abstraction for CYPs (Manchester et al. (1997) J. Am. Chem, 
Soc. 119, 5069-70). This abundant amount of data for both mammalian and bacterial 
isoforms strongly supports a common mechanism for multiple CYPs. 

A comparison of all the isotope effect values for the tested substrates except 
toluene revealed that the means were statistically different from each other at the 0.05 
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significance level. This can best be explained by the fact that the standard deviations 
reflect the precision of the method and do not take into account systematic and day to 
day errors. The design of intramolecular isotope effect experiments of substrates of 
the form R-CH2D specifically lends itself to very precise data. The magnitude of the 
5 intramolecular KIE is independent of substrate concentration, and no internal standard 
is necessary for the calculation of the results. 

In this regard, lEP's are enlightening. While 2B1 shows the smallest KIE for 
the P-OCH3 compound, it gives one of the higher means for other substrates. In fact, 
it is obvious from Figure 6 that no isoform consistently gives an lEP that is either high 

10 or low relative to other enzymes. 

Thus for the C YP en2ymes tested, the energetics of the active-oxygen species 
appear to be the same. 

In this example, all solvents were purchased from J. T. Baker, Inc. 
(Phillipsburg, NJ). Chemicals were purchased from Aldrich Chemical Co. 

15 (Milwaukee, WI) with the following exceptions: MTBSTFA (+1 % /er/- 

butyldimethylchiorosilane) was purchased from Regis Technologies, Inc. (Morton 
Grove, IL), and potassium phosphate dibasic and monobasic and magnesium sulfate 
were from EM Science (associate of E. Merck, Darmstadt, Germany), p- 
Chlorotoluene,/?-methylanisole, and /^•xylene were purified by distillation before use. 

20 /)-Bromotoluene was crystallized from ethanol and j:?-tolunitrile was crystallized from 
benzene:pet ether. Gibco BRL products (Gaithersburg, MD) were used for all cell 
culture methods. Biochemicals were purchased from Sigma (St. Louis, MO). p~ 
Xylem-CL^ a, a, a', a', a'-d6 (99.8% pure by GC) was purchased from Aldrich 
Chemical Company. 

25 Substituted Toluene-a-^Hj and -a-^Hj-a'-^Hi. Monodeuterated analogues of 

p-bromotoluene,p-chlorotoluene,p-methylanisole, and toluene were synthesized by 
the mesylation of the corresponding nondeuterated benzyl alcohols according to the 
same procedure described above for the trideuterated compounds, and /7-xylene-a- 
was synthesized with the same method from benzene dimethanol, except 
30 that the mole ratio of benzene dimethanol to the other reactants was 1:2. 

For p-tolunitrile-a-^Hi, the benzyl alcohol was first synthesized by the 
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reduction of /?-cyanobeiizaldehyde using sodium borohydride. For this procedure, the 
sodium borohydride (0.05 mol) was mixed with about 100 mL of ethanol and stirred 
under N2. The aldehyde (0.025 mol) was dissolved in a small amount of ethanol and 
added dropwise to the NaBKj solution. The reaction was monitored by TLC (silica 
5 gel: 50% ethyl acetate, 50% hexane). After about 1 .5 h, the mixture was filtered and 
the solvent was evaporated. The oil that remained was mixed with water and then 
extracted with CH2CI2. The organic layer was dried over MgS04 and evaporated to 
yield a white solid, which was recrystallized with hexane and assessed for purity by 
NMR and GC-MS. 

1 0 Preparation of HepGl-Expressed CYP, Individual C YP isoforms ( 1 A2, 2B 1 , 

2E1, and 2C9) were obtained using the HepG2 vaccinia expression system described 
previously (Gonzalez et al. (1991) Methods Enzymol 206, 85-92). Briefly, the 
recombinant vaccinia viruses containing cDNAs encoding a single CYP were 
propagated in himian TK~ 143 (human embryoblast) cells, which were then used to 

1 5 infect confluent flasks of hepatoma HepG2 cells (American Type Culture Collection, 
Rockville, MD). For each isoform, cells were scraped from 50 flasks, pooled, divided 
into 1-mL aliquots, and stored at -78 °C. Microsomal preparations from crude cell 
lysates were obtained immediately before use according to the following procedure: 
aliquots were thawed, mixed with 1 mL of 50 mM KPi buffer pH 7.4, sonicated with 

20 a Branson Sonifier Model 450, and centrifuged for 15 min at 300000g in a cold (4 ""C) 
Beckman TL-100 ultracentrifuge; pellets were reconstituted with 2 mL of KPi buffer 
pH 7.4, pooled, and homogenized in a cold room. The total protein concentration of 
the final enzyme-buffer mixture was typically 3-4 mg mL~* by Bradford assay 
(Bradford, M. M. (1976) Anal. Biochem. 72, 248-54) with bovine serum albumin as 

25 the standard. The typical low levels of P450 obtained using this expression system 
were undetectable by CO difference spectra by the method of Omura and Sato 
(Omura et al. (1964) J. Biol. Chem. 239, 2370-8). 

Cytochrome P450cam, Putidaredoxin, and Putidaredoxin reductase. 
Construction of expression plasmids for P450cam, Pd, and PdR as well as protein 

30 expression and purification have been described previously (Shimoji et al. (1998) 
Biochemistry 37, 8848-52). 
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Incubation Conditions and Isolation of Products. All incubations using 
HepG2-expressed CYP were run for 20 min in a shaking water bath at 37 °C and 
contained the following: 300-400 fig of protein suspension, depending on the 
isoform; 1 0 mM glucose-6-phosphate, 1 mM NADP'*', and 1 U glucose-6-phosphate 
5 dehydrogenase as an NADPH generating system; and 1400 U catalase in a total 
volume of 5 mL of 50 mM KPi buffer, pH 7.4, 2 mM MgC^. For P450cam, the 
assays were carried out at 37 °C for 20 min in a 1-mL reaction mixture in 50 mM KPi 
buffer, pH 7.4, 200 mM KCl, L5 ^M P450cam, 3 ^M Pd, 1.5 \xM PdR, and 0.3 mM 
NADH. Reactions were initiated by the addition of 50 \xL (HepG2-expressed CYPs) 

10 or 10 nL (P450cam) of a substrate solution in CH3CN, terminated upon addition of 3 
mL of CH2CI2, and 10 of internal standard in CH3CN was added, if applicable (see 
below). Samples were extracted three times with CH2CI2; the organic fractions were 
pooled, dried with anhydrous MgS04, and concentrated gently imder a stream of N2 to 
a volume of approximately 100 ^L, After reconstitution with 300 \xL (HepG2- 

15 expressed CYPs) or 100 jiL (P450cam) of CH3CN, samples were derivatized with 50 
\jLh of MTBSTFA + 1% /-BDMCS were heated overnight at 70 ''C. 

Enzyme incubations with monodeuterated toluenes and p-xylene-a-^Hi-a'-^Hi 
were performed in triplicate at a single substrate concentration of 0.25 mM. Controls 
for toluene-a-^Hi showed a small peak at the same retention time as product that 

20 originated from the derivatizing agent, which was quantified and subtracted from the 
observed amount of alcohol product. The kinetic parameters for metabolism of the 
substituted toluenes were determined by GC-MS analysis of the M-(dimethylsilyl) 
derivative of the product benzyl alcohols (see below). Control experiments using a 
mixture of deuterated and non-deuterated products showed that there was no isotope 

25 effect on product work-up. 

Product Analysis by GC-MS, All samples were analyzed with an HP (Palo 
Alto, CA) GC-MS 5890/5972 system, equipped with an HP 7673 A automatic 
injector. An HP-1 capillary column (25 m, 0.2-mm ID, 0.5-^m film thickness) was 
used for product analysis from all substrates except toluene and p-xylene, which were 

30 analyzed with an HPWax capillary colunm (30 m, 0.25-mm ID, 0.5-fim film 
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thickness). For the HP-1 column, the injector and detector temperatures were both set 
at 250 °C; for the HPWax column, the injector and detector temperatures were 230 
and 240 °C, respectively. All samples were injected in the splitless mode. Oven 
conditions for the HP-1 column were 50 °C for 0.5 min, 10 ^'C/min to 250 °C, and 2 
5 min at 250 °C except when /?-tolunitrile was used as the substrate, in which case a 
gradient of 1 5 °C/mm was used. For the HPWax column the oven conditions were 50 
°C for 0.5 min, 10 °C/min to 240 and 3 min at 240 °C. The MS was operated at 
an ionizing voltage of -70 eV and a 20-ms dwell time. Scans of authentic standards 
of the derivatized benzyl alcohols were used to find peak retention times and to verify 

10 products formed from enzyme incubations. Selected ion recording of the 

M-(dimethylsilyl) ion of products and internal standards was used to determine the 
amount of metabolism. The ratio of the protio and deuterio products was corrected for 
ion overlap and percent deuterium incorporation with Brauman's least-squares 
approach (Korzekwa et al. (1990) Biomed. & Environ. Mass Spectrom. 19, 21 1-7). 

15 Deuterium incorporation of substrates was measured with an HP GC-MS 5890/5970 
system with an HP-1 colimm using the same operating conditions as above except the 
MS ionizing voltage was adjusted to -10.0 eV. 

Data Analysis, Kinetic expressions were solved using the program 
Mathematica (Wolfram Research, Champain, IL). 

20 

Example 3 

A linear free-energy relationship (LFER) electronic model can describe CYP enzymes: 

Given the indication that the active oxygen is very similar in multiple CYP 
isoforms, a single electronic model can be used to describe (the corresponding) 

25 multiple CYP enzymes. One potential model is the use of linear free-energy 
relationships (LFERs), which are utilized in chemistry and enzymology as a 
predictive tool and as a probe of the electronic characteristics of the transition state of 
a reaction. The following tests a second hypothesis that LFERs can be used to predict 
relative rates of metabolism by CYP and to evaluate physical characteristics of CYP 

30 reactions. 
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Relatively few LFERs have been reported for CYPs. This is not surprising 
given the fact that the rate-limiting step in the catalytic cycle has been shown to occur 
prior to the step for substrate oxidation (White et al. (1980) Ann. Rev. of Biochem. 49, 
315-56; and Tyson et al. (1972) J, Biol Chem. 247, Slll-U). Good correlation of 
5 Fniax and feat values with either Hanunett a or oxidation-reduction potentials have 
been found for series of/?^7ra-substituted TViA^dimethylanilines (Burka et al. (1985) 
Am. Chem. Sac. 107, 2549-51; and Macdonald et al. (1989) Biochemistry 28, 2071-7) 
and para-substituted sulfides (Watanabe et al. (1980) Tetrahedron Lett. 21, 3685-8). 
These data have been used to support the mechanism of electron transfer from the 
10 lone pair of the heteroatom in the substrate to the active oxygen of CYP. Similarly, 
hydrophobic and electronic factors were found to be good descriptors of ^at values for 
a series of /^ara-substituted toluenes (White et al. (1986) Arch Biochem. Biophys. 
246, 19-32). 

LFERs correlate rates and equilibrium constants. The term free energy arises 
15 from the fact that equilibrium constants (K) and rates (k) are related to free energy 
(AG - -2.303RT logK and AG* = -2.303RT log(WkT)). LFERs can be used to 
observe the effects of substituents or solvents on the reactivity of a reaction. The 
degree and type of perturbation can provide mechanistic information about the 
chemical reaction. Rates for a reaction of interest are often compared to empirically 
20 determined constants. The most widely used empirical constants are Hammett 
constants. The Hammett constant a'^para is a physical parameter that indicates the 
relative degree of electronic perturbation that a substituent imparts in the para 
position. The Hammett constant CT"*^para is defmed by Eq. 1 1 where Kx is the acidity 
constant for a 

25 porj^-log^ (Eqll) 

para-substituted tert-cumyl chloride and Kh is the acidity constant for tert- 
cumyl chloride (in 90% acetone-H20). 

Given the above data (Example 2) supporting the concept that CYPs act with 
similar chemical mechanisms; the next question addressed was whether LFER can be 
30 used to predict the relative reaction rates of substituted toluenes. If relative rates of 



36 

BNSDOCID: <WO 000620SA1J_> 



wo 00/08205 



PCT/US99/17713 



product formation are masked by other steps in an enzymatic cycle, however, LFERs 
must be interpreted with caution. Kinetic information can be deduced, but not 
mechanistic information, if relative reaction rates are masked. Therefore, it is 
important to assess whether or not masking occurs with the substrates chosen for a 
5 LFER analysis. 

The first isoform chosen for kinetic studies was CYP2E1, since this is a low- 
Km isoform for small hydrophobic molecules (Guengerich et al. (1991) Chem. Res. 
Toxicol. 4, 168-79). The steady state kinetic constants Vmax and Km for hsnzylic 
hydroxylation were determined for five /7ar£i-substituted toluenes (p-X-Ph-CR^: X = 

10 OCH3, CH3, CI, Br, CN; Ph = phenyl) using expressed CYP2E1, which displayed 

typical Michaelis-Menton behavior (Table 6). The choices of para-substituents allow 
for a range of electronic perturbation at the benzylic group, thus a LFER with the 
Hammett constant d^p was originally explored (Hansch et al. (1979) Substituent 
Constants for Correlation Analysis in Chemistry and Biology, John Wiley & Sons, 

1 5 New York). Hammett constants have been shown in other cases to correlate well with 
rates of hydrogen atom abstraction from substituted toluenes (Kim et al. (1985) 
Tetrahedron Lett 26, 891-4; Sakurai et al. (1967) J. Am. Chem. Soc. 89, 458-60; and 
Heberger, K. (1994) J. Phys, Org. Chem, 7, 244-50). 
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Table 6. Kinetic Constants for the Formation of Para-Substituted Benzyl 
Alcohols from the Corresponding Para-Substituted Toluenes Using Expressed 



CYP2E1 



para 
substituent 




S.E." K„ax 




S.E. V/K 




n^ 


CN 


61.7 


0.4 


0.014 


0.001 


1 


6 


Br 


58,5 


1.3 


0.043 


0.005 


0.979 


12 


CI 


50.4 


0.8 


0.033 


0.002 


0.998 


6 


CH3 


58.6 


1.2 


0.022 


0.002 


0.990 


11 


OCH3 


41.0 


1.1 


0.017 


0.002 


0.983 


11 



Br-d3 


54.8 


0.8 


0.027 


0.002 


0.997 


6 


Cl-d3 


31.4 


1.0 


0.020 


0.003 


0.991 


6 


*CH3-d6 


55.7 


0.4 


0.019 


0.001 


1 


6 


OCH3-d3 


18.0 


0.3 


0.007 


0.001 


0.996 


6 



5 ° Kmax aild V/KvfQxe determined by nonlinear regression. * Units of nM min*. 

^Standard error of the indicated constant. Units of min~\ * Correlation coefficient 
for the nonlinear regression fit to kinetic data. -^Number of points used for the 
determination of kinetic constants J^max and V/K. Experiments were performed with 
nondeuterated or trideuterated (at the benzylic position) substrates, except for 
1 0 /?-xylene. ^/7-Xylene-a, a, a, a', a', a '-de was used as the deuterated analogue of 
77-xylene. 

Only poor correlation was found in a plot of Fmax versus a^p (Log Kmax = 
0.1 1 ar% + 1.73; r^ = 0.588) as well as other physical constants (data not shown), and 
1 5 poor correlation with V/K and various parameters. In addition to the poor correlation 
of Kmax and a'^p, the slope is positive, which is unexpected based on the rates of H- 
atom abstraction from toluenes using various oxygen radicals, which yield negative 
slopes (Kim et al., Ibid\ and Heberger et al., Ibid,\ A positive slope indicates that a 
transition state is stabilized by electron withdrawing substituents. This would 
20 normally be interpreted to mean that the transition-state is negatively charged. 

Three possibilities were considered for the explanation of the poor 
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correlations: 1) electronic parameters do not determine CYP2E1 substrate reactivity, 
2) the kinetic constants (Fmax and V/K) are masked or partially masked, and 3) the 
transition state is negatively charged. The first two possibilities were analyzed with 
the appropriate methodology and the results are discussed in detail below. 
5 Correlation of electronic parameters of substrates with KDIE, First, the 

question of whether or not electronic factors of CYP substrates have any observable 
effect on the oxidation step was addressed. The natural log of the intramolecular 
KIEs for the six toluenes listed in Table 5 above were plotted as a function of a% for 
each /?-substituent. The substituent and corresponding ct'*^p value are (respectively): 

10 OCH3, -0.78; CH3, 0.31; H, 0; CI, 0.11; Br, 0.15; CN, 0.66 (Hansch et al., Md), In 
this series of substituted toluenes, it is proposed that the force constant between the 
benzylic carbon and the abstractable hydrogen is different for each substituent, and 
varies as a function of a'^p. The data in the poor correlation plot of Finax versus a'^p 
(described above) support this idea. A KIE plot (of the six toluenes listed in Table 5 

1 5 above versus a'^p for each /7-substituent) demonstrated values that lie on the ascending 
slope of the parabolic curve described above in Example 2. It is concluded that the 
transition-state is sensitive to electronic effects, and that electron-donating 
substituents stabilize the transition-state. 

KIEs and Linear Free-Energy Relationships. When LFERs are sought, 

20 experiments should be performed that definitively answer the question: is the rate of 
product formation masked (is product formation limited) by a rate determining step in 
the catalytic cycle that is not the chemical step of interest? Fortunately, innovative 
methods for investigation of masking in enzymatic systems were developed by 
Cleland and Northrop (Northrop (1975); and Cleland, W, W. (1975) Biochemistry 14, 

25 3220-4). With this approach, KIEs are measured on the enzymatic parameters V^ax 
and V/K, Isotope effects on Kmax and V/K are expressions for the isotope effects on 
rates of product formation with the limiting conditions of high and low substrate 
concentrations, respectively. 

Reported herein are deuterium isotope effects on Vmax and PTST which are 

30 designated ^Kand ^(V/K). The expressions for both '^Kand ^(V/K) contain the rate 
constant for the isotopically sensitive step, but they differ in the following ways: only 
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contains the rate constant for dissociation of product from enzyme, and only 
\V/K) contains the rate constants for binding of free substrate and enzyme, and 
dissociation of the first ES complex to free substrate and enzyme. The relative 
magnitudes of these steps as well as others to the isotopically sensitive step determine 

5 the degree to which the intrinsic KIE (the KIE associated solely with the bond 

breaking step, /.e., kH/ko) will be expressed in ^Kor ^(V/K), Therefore, comparisons 
of ^Fand ^{V/K) isotope effects to an intrinsic KIE can provide information regarding 
which steps, if any, are involved in masking the relative rates of product formation. 
These types of studies can reveal if LFERs can be used to make mechanistic 

10 conclusions. It should be noted that even if a linear relationship is seen between In 

rates and descriptors for a series of compounds, it does not mean that the relative rates 
are unmasked. 

Intermolecular KIEs were used to investigated whether or not the poor 
correlation of Vmax and ct% was caused by masking of K^ax by some other (slower) 

15 step in the catalytic cycle. Most enzyme catalyzed reactions contain multiple steps, 
and it is not uncommon for masking of reaction rates to occur. For CYP reactions, it 
is generally accepted that the rate limiting step is transfer of the second electron from 
P450 reductase to P450, a step that occurs prior to formation of product. 
Consequently, experimentally measured kinetic constants may be partially masked 

20 due to the presence of the slow electron transfer step unless branching to an alternate 
position is fast enough to immask the intrinsic value for the kinetic constant tmder 
investigation. In summary, information gained from an enzymatic rate me£isurement 
is mechanistically meaningless if the rate is not that for the reaction step(s) of interest 
To ascertain if the Kmax values from Table 6 were masked, intermolecular 

25 isotope effects were measured and compared to the corresponding intramolecular KIE 
from Table 5. Isotope effect experiments of the noncompetitive intermolecular design 
were performed in order to evaluate the two kinetic constants Fmax and V/K, Table 7 
shows that the observed KIE on both Kmax and V/K for four of the /^-substituted 
toluenes listed in Tables 5 and 6 range from 1 to 3, while the intramolecular KIE for 

30 the corresponding toluenes in 2E1 range from 4 to 8. Since the intermolecular KIEs 
are smaller than the intramolecular KIEs, it can be concluded that both ^Fand ^(V/K) 
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are masked. The fact that ^Vpi = ^(V/K) pi values are ahnost identical support the fact 
that the external conunitments to catalysis are small (see Scheme 3 and corresponding 
discussion). 



Table 7. Comparison of Isotope Effects on Benzyl Alcohol Product Formation 
for Various Para-Substituted Toluenes Among Three Different Types of Isotope 
Effect Experiments 



Intramolecular' 



para 
CYP substituent 



kH/kD 



Intermolecular 
Noncompetitive ^ 



Intermolecular 
Competitive* 



2E1 



6.75 ± 0.03 1,07 ± 0.04 1.58 ± 0.01 



2.06 ± 0.07 



CH3'' 5.45 ±0.05 1.05 ±0.03 1.19 ±0.01 



1.18 ±0.03 



OCH3'' 4.24 ±0.01 2.28 ±0.10 2.60 ±0.01 2.57 ±0.02 



cr 



6.75 ±0.02 1.61 ±0.01 1.62 ±0.01 



3.18 ±0.62 



nd 
nd 



3.21 ± 0.02 
2.47 ±0.10 



2B1 cr 8.06 ±0.12 nd^ 

1A2 cr 7.06 ±0.02 nd 

" The numbers after the ± sign are propagated errors. * Intermolecular 
competitive values are expressed as mean ± standard deviation of three samples. The 
10 substrate analogues used for experiments of intermolecular design were either ''do and 
da, or ''do and de; only benzylic position(s) were deuterated. ^Not determined. 

Intermolecular isotope effects, while usually not providing significant 
mechanistic information since the intrinsic isotope effect is not observed, can provide 
15 information about the kinetic mechanism. For example, this type of experiment is 
traditionally used to determine if a bond-breaking step is rate-limdting. An 
intermolecular noncompetitive KIE experiment involves the independent 
determination of Fmax and V/K using protonated and isotopically labeled substrates, 
for example, R1-CH3 and R1-CD3. For the calculation of a KIE, the kinetic constant 
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(^max or V/K) pertaining to unlabeled substrate is divided by the corresponding 
constant for labeled substrate. 

Scheme 2 shows a potential kinetic mechanism that describes an 
intermolecuiar noncompetitive KIE, and equations 12 and 13 are the expressions for 
5 ^Kand "^{V/K) obtained for this scheme. Substrate (Sh or Sd) combines with enzyme 
(E) reversibly to form the corresponding ES complex (ESh or ESd). (Subscript 'H' or 
*D' signifies the involvement of a compound or complex in the pathway for 
abstraction of hydrogen or deuterium, respectively). Formation of an * activated' 

10 Scheme 2 

E + Sh [EShI [ESHlact-^ [EPhI E + Ph 
k2 

E+Sd [ESd] [ESEOact-^^ [EPd] E + P^ 
^2 



'-SH 



''5H 



15 



1+ 



(Eq. 12) 



^{V/K) = \ (Eq.l3) 

20 complex ESh* or ESd*, from ESh or ESd, respectively, is depicted here as an 

irreversible step that pccurs prior to the isotopically sensitive step. This irreversible 
step is consistent with the understanding of the kinetic mechanism of C YP, It is a 
kinetic representation of a collection of steps that include: the first electron reduction 
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10 



15 



20 



(Fe -> Fe oxygen binding and reduction; the second electron reduction of 
oxygen; and the heterolytic cleavage of the peroxy anion, a postulated step which is 
presumed to be irreversible. Product is formed at the isotopically sensitive step, 
producing an enzyme-product complex (EPh or EPd), v^hich is followed by 
dissociation of the corresponding product from the enzyme. Rate constants k}, fe, 
and ky may each actually represent more than one step, and inclusion of additional 
steps increases the complexity of the KIE equations but does not modify this analysis. 
When there is no pathway for alternate product formation, which is typical for most 
enzymes, the value of '^F depends on the relative magnitudes of Jta and kj to ksH- As 
the rate constants fe and fe decrease relative to ksn, becomes masked (see Eq. 12). 
^(V/K) for the noncompetitive KIE with the kinetic mechanism shown in Scheme 2, 
which has no branching pathway from ESh* or the corresponding ESd* position, is 
completely masked (Eq. 13). 

While the ^(V/K) is completely masked given Scheme 2, a branched pathway 
has the potential to unmask both ^Kand "^(V/K), This is illustrated using Scheme 3, 
which represents the formation of two different products from one substrate, either 
deuterated or non-deuterated. Metabolism at an alternate, isotopically insensitive 
step, fe, is in competition with the isotopicaUy sensitive step, fe. Rate constants ki 
and fe represent substrate binding to and debinding from the enzyme. Rate constant 
fe is shown as an irreversible step that represent steps that occur after substrate 
binding, up to and including oxygen activation. These steps lead to formation of 
'activated complexes' (ESI* and ES2*) between the enzyme and the substrate in the 
relative positions 1 and 2. The active oxygen species is stable enough to permit 
translational and conformational changes of substrate in the active site of the enzyme. 



Scheme 3 



kl 



E + S 



k2 



kl 




ES2* 



EP2 



E + P2 
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_ ^7D 



^7H^6 



^7H 



2 I ^7H^6 



(Eq. 14) 



^7H ^ 
^7D 




_^ ^7H 
^5 . 




2^3 ^7H J ^3^7H J ^3^7H | '^^:3^S^1H 
_ ^13^7D ^1 1^5 ^13^5 ^1 1^5^ 


1 + 


^7H^6 
_ ^5^9 


, ^7H 

^5. 


+ 


2^3 ^ A:3^7H 1 ^3^7H j 2^3^6^7H 
_^13 ^U^5 ^13^5 ^11^5^9 _ 





(Eq. 



5 15) 



Consequently, it is possible for an equilibrium to be established between ESI * 
and ES2* (see ks and ke). This is the key to unmasking KIEs. When deuterium is 
substituted for hydrogen at the site of metabolism, the rate constant for product 

1 0 formation is reduced. This can be expressed as A:7d < km, where ^70 represents the 
rate constant for oxidation at a deuterated position on a substrate, and km represents 
the rate constant for oxidation at the same position on a nondeuterated analogue of the 
substrate. Deuterium substitution will cause the concentration of ESI d* to increase 
relative to the concentration of ESI h* unless ESId* is in equilibrium with another 

15 complex, for example, the ES2d*. If an alternate pathway does not exist, then 
[ESId*] > [ESIh*], and an isotope effect is masked. 

Oxidation at an alternate, imlabelled position on a substrate with the formation 
of a second product (P2) can unmask a KIE. This is accomplished by the redirection 
of "excess" ESId* towards ES2d* which allows [ESId*] = [ESIh*], which is 

20 necessary for the expression of a maximum, or unmasked, isotope effect. According 
to convention, a substrate with a low forward internal commitment to catalysis will 
branch from ESI towards ES2 (Scheme 3), and isotope effects on "^Kand ^{V/K) will 
be unmasked, depending on the magnitude of the conamitment terms. The extent to 
which ^{V/K) is unmasked depends on the relative values of ks and k9 to kvH (Eq. 14). 

25 As the rate of substrate rotation and branching increase relative to the rate of PI 
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formation, the rate of unmasking increases. The expression for a KIE on Kfor 
Scheme 3 is quite complex (not shown) but can be simplified to Eq. 15 with the 
following two assumptions. If it is assumed that reduction of oxygen is rate limiting 
(see prior disciassion), and is significantly smaller than: (i) rotation of substrate in the 
5 active site (^3 « ^5, fe); and (ii) the rates of product formation k^H and k^, the 

expression for ^Kpi is simplified to Eq. 15. If the external reverse commitment to 
catalysis is small (kn, kis » ^3), i.e., product release is not rate limiting, and Eq. 15 
reduces to Eq. 14 (^Fn = ^{F/K)pi). 

CYP has the potential to function as an oxidase, whereby molecular oxygen 

1 0 undergoes a four electron reduction, with the overall formation of two molecules of 
water. Formation of the second molecule of water acts as a branch point leading from 
the ES complex in Scheme 3, and precludes substrate oxidation. As a consequence, if 
branching to water is fast relative to product formation, KIEs (or rates of product 
formation) can become unmasked- 

1 5 The rate equations 12 and 13 for an intermolecular noncompetitive KIE 

without branching to alternate product (see Scheme 2) can be used to interpret these 
results. As the relative ratio of ^sh to ^3 and/or k^ increases, the degree of masking of 

increases. An irreversible step preceding the isotopically sensitive step, such as 
formation of the active oxygen species, can mask '^Kif it is slow enough. In addition, 

20 an irreversible step following product formation, such as product release, can produce 
the same effect, depending on its magnitude. For a noncompetitive scheme with no 
branching from the activated complex ESh* or ESp*, ^(V/K) will be equal to 1 . Since 
the ^{V/K) KIEs observed are larger than one, branching to an altemate pathway (not 
shown in Scheme 2) may account for the partially immasked values of 1 .2 - 2.7 

25 observed (see Eq. 13 and relevant discussion). 

Branching to an altemate oxidation site(s) on the substrate, or through the 
pathway leading to water formation have been shown to immask isotope effects for 
various compounds. With most of the substrates used in this study, it is possible that 
branching to water serves to slightly unmask ^(V/K). Benzylic hydroxy lation 

30 accounts for 70% or more of the total oxidation of toluene, depending on the source of 
CYP, and -95% of /7-xylene metabolism. Ring hydroxylation is energetically less 
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favorable for p-bromotoluene, p-chlorotoluene and p-tolunitrile, as compared to 
benzylic hydroxylation. In these cases, since metabolism at the alternate position on 
the substrate is slower than metabolism at the benzylic group, it cannot fully unmask 
an isotope effect. The data in Table 7 indicate that at least for p-methylanisole and p- 
5 chlorotoluene the intermolecular isotope effects at the benzylic position are slightly 
unmasked. One likely interpretation is that a small amount of branching to water is 
responsible for this unmasking- Furthermore, since ^(V/K) and isotope effects are 
nearly identical, the masking must be due to high internal commitments to catalysis. 
The distinction between internal and external commitments to catalysis are described 

10 by Northrop (Methods Enzymol 87, 607-25, 1982). 

In order to confirm the intermolecular noncompetitive KIE data, 
intermolecular competitive \V/K) were measured for the same substrates using 
CYP2E1 (Table 7). Isotope effects of the intermolecular competitive design are 
simpler to perform as compared to experiments of noncompetitive design since only 

15 one concentration of each substrate (deuterated and nondeuterated) is used. In this 
type of experiment, both substrates are added to the same incubation mixture at a 
single (usually equal) concentration, and a competition between the deuterated and 
nondeuterated compounds is established. Only isotope effects on V/K are measured 
with experiments of this design. 

20 The results of the competitively determined ^{V/K) for four substituted 

toluenes are listed in Table 7. The values of 1.2 - 3.2 indicate that the KIE are 
masked, when compared to the intramolecular KIE values for the same substrates, 
which range from 4 to 7. For /?-xylene and /?-methylanisole, the ^(V/K) KIE 
determined competitively and noncompetitively are the same. The KIEs for p- 

25 bromotoluene determined by competitive and noncompetitive experiments are 1.5 and 
2,1, respectively, and for p-chlorotoluene the values are 1 .6 and 3.2. The 
discrepancies in ^(V/K) for /7-bromotoluene and /7-chlorotoluene can be interpreted in 
terms of the kinetic constants that define the schemes for the two types of 
experiments. In experiments of the competitive design only, steps that pertain to 

30 debinding of both deuterated and nondeuterated substrates from one (or more) of the 
ES complexes can potentially urunask (or at least partially immask) a KIE, depending 
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on the magnitudes of the rate constants for these steps relative to others. The reason 
for the substrate dependence on the differences in ^(V/K) for the toluenes is not 
known at this time, but it is certain that KIEs are masked in both types of experiments. 
Competitive KIE experiments were performed for p-chlorotoluene using 
5 expressed 2B 1 and 1 A2 to address whether or not masking of ^(V/K) occurs in other 
CYPs (Table 7). The K/ii: isotope effects were 3,1 (2B1) and 2.5 (1A2), as compared 
tointramolecular KJEof 8.1 (2Bl)and7.1 (1A2). The results imply that the same 
degree of masking occurs in all three isoforms (for /?-chlorotoluene). 

A number of LFER with CYP can be found in the literature. Burka et al. 

10 (Proc. Natl Acad. Set USA 80, 6680-4, 1983) found a linear correlation between V^ax 
for aromatic hydroxy lation and a"*", and argued for an initial electron transfer 
mechanism directly from the benzene ring, for a small nimiber of substituted 
monohalobenzenes. Rates of CYP-mediated A^-dealkylation was correlated with 
oxidation-reduction potential, in accordance with Marcus theory, according to 

15 Guengerich, et al. (J, Biol Chem. 271, 27321-9, 1996). Rates of S'-oxygenation by 
CYPs were correlated with Hammett to give a negative slope, which the authors 
used to argue the development positive charge in the transition state (Watanabe et al. 
(1980) Tetrahedron Lett. 21, 3685-8). White and McCarthy found a LFER with a 
series of substituted toluenes {Arch, Biochem. Biophys. 246, 19-32, 1986). Both the 

20 electronic parameter Hammett sigma, and a lipophilicity parameter correlated well 

with their data. However, isotope effects for a of substituted A^,A^-dimethylanilines by 
Dinnocenzo, et al {J. Am. Chem. Soc. 115, 71 11-6, 1993) showed that isotope effects 
were masked, and therefore rates, were masked. The LFER and KIE data also 
supports the theory that rates and isotope effects are for CYP reactions are masked. In 

25 addition, the isotope effect for toluene that was measured by White and McCarthy 
was 2, which indicates that the rates were masked {Arch. Biochem. Biophys. 246, 19- 
32, 1986). In conclusion, it has been shown that the chemical step is masked for 
benzylic hydroxy lation of a series of substituted toluenes. This is likely to apply to all 
LFERs that have been done on CYP to date. 

30 In this example, chemicals and reagents were obtained as described in 

Example 2. 
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Synthesis of Compounds. Para-Substituted benzyl-a-^ H2 Alcohols , Methyl 4- 
bromobenzoate, methyl 4-chlorobenzoate, methyl 4-cyanobenzoate, methyl 4- 
methylbenzoate, and methyl 4-methoxybenzoate were used as starting materials for 
the synthesis of the corresponding para-substituted benzyl alcohols according to the 
5 following procedure. Lithium aluminum deuteride (1 .2 mol) was suspended in ether 
and cooled to 0 °C in an ice bath while under N2. The substituted benzoate (1 mol) 
was dissolved in a small amount of ether and added drop-wise to the LAD suspension 
while keeping the temperature at 0 After the addition, the mixture was stirred at 0 
°C for an additional hour, or until complete by TLC (silica gel: 80% hexane, 20% 

10 ethyl acetate). The xmreacted LAD was decomposed, and the mixture was filtered, 
dried over MgS04, and evaporated to dryness. The crude product was purified by 
recrystallization and/or distillation, and purity was assessed by NMR and GC-MS. 
All NMR and GC-MS data was consistent with commercially available compounds. 
Para-Substituted Toluenes-a, a, a-^/fj. p-Bromotoluene,/7-chlorotoluene, 

15 and/>-methylanisole were synthesized firom the corresponding j>-substituted benzyl-^/? 
alcohols according to the following method. Methane sidfonyl chloride (0.1 1 mol) 
was dissolved in a small amount of CH2CI2 and added drop-wise to a solution of the 
^substituted benzyl-di alcohol (0.10 mol) and triethylamine (0.15 mol) in CH2CI2 
that was cooled to -78 °C in a dry ice-acetone bath. The reaction was stirred for 

20 another 30 minutes at -78 °C. Because the mesylates are generally unstable, they 
were not isolated, and the procedure was continued by adding ice water, separating 
the layers, and extracting the aqueous layer with CH2CI2. The organic layers were 
combined, washed with saturated sodium bicarbonate, and dried over MgS04, and the 
solvent was removed under reduced pressure. The oil that remained was reconstituted 

25 in ether and reduced with LAD using the same procedure as described for the 

reduction of the benzoates. The products were purified by chromatography (silica 
gel: 100% CH2CI2), distilled under vacuum using a Kugelrohr, and, when necessary, 
purified again with a neutral alumina column and 100% hexane to remove any 
residual benzyl alcohol. Any residual benzyl alcohol was detected by select ion 

30 monitoring using GC-MS as described below. 

Preparation of HepGl-Expressed CYP. See Example 2. 
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Incubation Conditions and Isolation of Products. As in Example 2 and the 
following. Upon termination of incubations by the addition of 3 mL of CH2Cl2,l 0 
nmol of the appropriate internal standard (in 10 fiL of CH3CN) was added, if 
applicable (see below). 
5 Determinations of F^ax and Km for Do or D3 substrates were made using at 

least six samples, with six different substrate concentrations in a range between 0.5 
and 100 |iM, and 10 nmol of the appropriate para-substituted benzyl alcohol as the 
internal standard. Product formation from Do substrates was measured with the 
corresponding para-substituted benzyl alcohol dideuterated at the benzylic carbon as 
10 the internal standard; for D3 substrates the internal standard was the corresponding 
nondeuterated benzyl alcohol. Competitive experiments were performed in triplicate 
with mixtures of Do and D3 substrates, each at a final concentration of 0.25 mM. 
Data Analysis. See Example 2. 

15 Example 4 

Rapid substrate rotation in CYP mediated hydrogen atom abstraction reactions: 

While the regioselectivity of an uncatalyzed reaction reflects the energetic 
differences for oxidation of the various positions in a substrate molecule, an enzyme 
catalyzed reaction the substrate molecule can be bound in a specific fashion to the 

20 enzyme, and bias the regioselectivity of reaction to favor a higher energy pathway 

than the uncatalyzed reaction. The separation of the reactivity differences intrinsic to 
the substrate from the steric interactions imposed by the enzyme is a challenging 
problem. This problem becomes particularly important when developing 
computational models for an enzymatic system, such as CYP, since the intrinsic 

25 reactivity differences need to be addressed with quantum chemistry, while the steric 
factors need to be addressed with molecular mechanics/dynamics. In this example, a 
method for determining if the intrinsic reactivity differences are free of CYP steric 
interactions is described using a combination of isotope effects and the 
regioselectivity of the reaction. 

30 A set of small substrates were used for the development and testing of 
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methods to determine factors that govern regioselectivity. The substrates o- andp- 
methylanisole and a-chloromethyl-p-xylene (Figure 7) are relatively small, and may 
be free from steric constraints imposed by the enzyme active site for most 
CYPisofr>rms. Upon comparing regioselectivity and KIE data among multiple 
5 isoforms, a rank order of reactivity for the fimctional groups found in the substrates 
shown in Figure 7 was established. 

Thus a process has been established by which computational and experimental 
data with different energy scales can be combined to give a unified model of 
reactivity. This approach allows the prediction of regioselectivity and rates of drug 

10 metabolism reactions. The computational methods that can be developed from this 
technique allow the extension of established computational methods for hydrogen 
atom abstractions to aromatic oxidation and other mechanistically distinct CYP 
mediated reactions. Furthermore, this method will allow the incorporation of 
computational and experimental models for other reactions involved in drug 

15 metabolism. This approach will assist in the assessment of metabolic properties of 
new drugs and in the redesign of promising compounds Avith better metabolic 
characteristics. 

The schemes and the corresponding arguments in this example provide the 
fundamental concepts for the formulated method which allows the separation of 

20 electronic factors from steric factors for CYP-mediated enzymatic reactions. Figure 
8 illustrates the CYP-mediated formation of two potential products from one 
substrate. Oxidation at the fimctional group near the tail of an arrow leads to the 
formation of the product near the head of the corresponding arrow. The goal is to 
answer the following question: does the observed ratio of products reflect only the 

25 intrinsic electronic differences between the two functional groups, or does the protein 
interact with or constrain the substrate, such that the product ratio reflects both 
electronic and enzyme-based steric factors? 

Scheme 4 illustrates a kinetic model for the formation of two products from 
one substrate. The details of each step and the nomenclature have been described in 

30 Example 3. It has been inferred from isotope effect experiments using various 

substrates and numerous CYP isoforms that the lifetime of the active oxygen permits 
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some degree of substrate reorientation or movement after formation of the active 
oxygen species. Thus, ESI * and ES2* can interchange if fe, and > zero. In 
addition, an equilibrium can exist between complexes ESI * and ES2*, which is 
defined by the equation {kslk^l{ks'lh'). 



Scheme 4 



V 

ki k3 
E + S ^ ^ ES ^ ES* 



k7 kii 
ESI* EPl ^ E + PI 



k2 

k5'^ ES2* > EP2 E + P2 

M = ^iZmil (Eq. 16) 



10 



Pill 



obs 



1 + 



Jk6 



(Eq. 17) 



Equation 16 shows the relationship between the observed product ratio 
(([FI]/[P2])obscivcd) and the intrinsic ratio for the relative rates of product formation, 
k7/kg. {The intrinsic ratio kj/kg can be likened to the ratio that would be observed if a 

15 chemical reaction occurred between a hypothetical [FeO]^^ species and the substrate, 
in the absence of protein.) The value for the observed product ratio depends on the 
intrinsic ratio for the relative rates of product formation and [ESI ♦]/[ES2*]. Equation 
16 can be rearranged to Eq. 17 upon solving for the steady-state concentrations of 
ESI * and ES2* {and assuming that Keg = 7). In order for the observed product ratio 

20 to accurately reflect the intrinsic differences in rates of metabolism, two criteria must 
be met. Criterion 1) The ESI * and ES2* complexes must be able to interchange 
rapidly relative to product formation, that is, k^» kTmEq. 17. Criterion 2) The 
equilibrium constant Kcq for interchange of ESI * and ES2* must be equal to 1 , that is. 
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[ES1*]/[ES2*] = 1, or ks = ^5 and h' = fe'- 

Criterion 1 can be evaluated by assessing rates of rotation with KIEs. First, 
intcrmolccular competitive isotope effects were measured using sets of substrates that 
have one site of metabolism that is isotopically substituted and has at least one other 
5 site of metabolism. Scheme 5 shows a mechanism that describes an intermolecular 
competitive KIE for substrates that have two oxidation sites. The details of each step 
and the nomenclature have been described in Exampl 3- In Scheme 5, SI represents 
nondeuterated substrate and S2 represents deuterated substrate. Oxidation occurs 
irreversibly at either the isotopically labeled position (through steps A:7h for ESIh*, 
1 0 and kjo for ES2d*), or the alternate, unlabeled site (through steps k\ 1 for ES 1 alt* and 
Jtii for ES2alt*» assuming the chemistry of this step is not isotopically sensitive). 



15 



Scheme 5 



ESUlt*— ^ EPIalt — ^ PIalt + E 



ESI ^ > ESI* ESlH* EPI -^5-^ PI + E 

^1 ^ ks 




k2 



E + SI + S2 



ES2 ^ > ES2* ^=F=^ ES2d* ^^> - EP2 > P2 + E 

ks 

ku kn 
ES2alt* — EP2alt P2alt + E 



kiH km 

«'>^= , km ^^-^^^ 
1 + — 

^6 



JP1ajlt]/[P2alt]J 



Equation 18 is the kinetic expression that describes the observed, combined 
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intermolecular competitive KIE for Scheme 5 when fe' > 0, and assuming Kcq — 1 for 
both pairs of ESI* and ES2* complexes. The step fe' is a step that leads to alternate 
ES complex formation. The combined intermolecular competitive KIE is the ratio of 
the KIE for product formation at the isotopically sensitive position ([P1]/[P2]) to the 
5 KIE for the formation of the altemate, isotopically insensitive product 
([PI alt]/[P2alt])- Equation 18 shows that when ke » ^th, the combined 
intermolecular competitive KIE is equal to the intrinsic KIE (^7h/^7d) for that 
compound. 

Equation 1 8 also describes how altemate product formation can unmask 

10 isotope effects and relative reaction rates for CYP reactions. The formation of a 
second molecule of water can also serve as an altemate product and therefore can 
unmask isotope effects. If water formation is increased upon deuteration of substrate, 
a large intermolecular isotope effect could be observed even if exchange is not rapid 
between the two ES* orientations. Rapid exchange can be confirmed by a 

1 5 comparison of the overall product formation from deuterated and protio substrates, 
independently. The degree of unmasking can be made by comparing the value for the 
intramolecular and intermolecular KIEs (described below). 

Next, intramolecular isotope effects, which can serve as very close estimates 
of intrinsic isotope effects (^h/Ad) when substrates are chosen that contain at least two 

20 symmetric sites, with at least one site (but not all sites) isotopically labeled, were 
measured. An example for the formula for this type substrate is R-CH2D. A 
comparison of an intramolecular isotope effect to the combined iptermolccular 
isotope effect allows the assessment of the magnitude of kjn/ke in Eq. 18. If the value 
for the intermolecular KIE = the value for the intramolecular KIE, and total product 

25 formation does not decrease upon deuteration of substrate, it can be concluded that h 
» k7H, and that the substrate undergoes rapid rotation in the enzyme active site. If 
the combined isotope effect and the intrinsic isotope effect are not equal, it is possible 
that steric interactions hinder rotation and influence regioselectivity. 

To evaluate criterion 2, the Kcq of multiple CYP isoforms are compared. The 

30 assessment, of the regioselectivity results, is whether or not the equilibrium constant, 
Kcq, for the two "activated" complexes ([ESI*] and [ES2*] in Scheme 4) is equal to 1. 
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It is possible that a preferred binding mode, imposed by the enzyme, exists for one of 
the two (or more) orientations of the substrate in the active site. If there was a 
preferred binding mode, then the concentrations of the activated complexes would not 
be eqiial. If [ESI *]/[ES2*] in Eq, 17 is not 1, then it is not possible for the observed 
5 regioselectivity to equal the intrinsic differences in rates of metabolism at the two 

positions (ki/kg). This concept would appear to be in violation of the Curtin-Hammett 
Principle if equilibrium between these two complexes were rapid. The Curtin- 
Hammett Principle states that if the ground state conformations are in rapid 
equilibrium only differences m the transition state energies of the two reaction 

10 pathways, not groimd state energies, will dictate which product is formed more 

rapidly. However, it is assumed that poor binding in the ground state would result in 
poor binding in the transition-state, and would introduce a steric effect on the 
resulting product ratios. Thus, only enzyme-substrate complements with K^q =^ 1 will 
provide good reactivity data. An indirect method to assess if K«q = 1 is to compare a 

1 5 number of individual isoforms to test whether regioselectivity is constant. 

In summary, if a combined intermolecular competitive KIE for a substrate 
that has two metabolizable positions is equal to the intramolecular KIE for that 
compound, then the two products are formed from kinetically indistinguishable 
complexes as long as water formation is not responsible for immasking the 

20 intermolecular isotope effect. If regioselectivity is the same for multiple isoforms, 
there is strong evidence that K«, = 1. Hence, regioselectivity results for a large 
number of isofonns provides a measure of the intrinsic differences in reactivity for the 
two functional groups as long as K^q = 1 for the kinetically indistinguishable 
"activated" complexes, and water formation isn't faster than product formation. If an 

25 isoform shows deviant regioselectivity results compared to a large number of 

isoforms, it is likely that enzyme-substrate interactions affect the regioselectivity for 
the deviant isoform, and that K^q ^ 1 . 

Regioselectivity of Ortho- and Para-Methylanisole. Benzylic hydroxylation 
(BzH) of o- and p-methylanisole (o- and p-MA) by C YP yields o- and p- 

30 methoxybenzyl alcohol, respectively. O-Dealkylation (ODAlk) of o- and ;7-MA 

yields o- and />-cresol, respectively (Figure 9). Metabolism of the ortho isomer also 
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yields 4-inethoxy-3-methylphenol (4M3MPh) from aromatic hydroxy lation. 
Regioselectivity was measured for dO-p-MA using expressed CYPs 1 A2, 2B1, 2B6, 
2C8, 2C9, 2EI and 4B1, and the results are shown in Table 8 as the ratio of CYP- 
mediated BzH to ODAlk. For all isoforms that were tested, BzH is favored over 
5 ODAIk. The range of regioselectivity is from 1 .3 to 20, which corresponds to a free 
energy of activation range (AAG*) of 0.5-1.8 kcal/mol. 
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The range of regioselectivity for o-MA, when expressed as BzH:ODAlk was 1 .2- 
3 for CYPs 1A2, 2B1, 2C9, and 2E1 (AAG* 0.1-0.65), and for 2B6 and 4B1 the 
regioselectivities were 0.52 and 0.7, respectively (Table 9). The regioselectivity range, 
when expressed as the ratio of BzH to aromatic hydroxylation (Arom) was 0.3-0.6 (AAG^ 
5 -0.7 to -0.3) for isoforms 1A2, 2B1, and 2E1, and 1.6-22 2E1 (AAG* 0.3-2) for isoforms 
2B6, 2C9, and 4B1. Isoform 2C8 produced about the same amount of product from all 
three postitions. Ring hydroxylation products accoxmted for approximately 30-60% of 
total products for all isoforms except 2B6 and 4B1, which gave 15 and 2% phenol 
product, respectively (Table 9). 

10 

Table 9. Metabolite Profile for o-Methylanisole 
Expressed as Mol-% 



CYP 


o-MBA 


o-CR 


Phenol 


1A2 


34 


12 


54 


2B1 


19 


18 


63 


2B6 


29 


56 


15 


2C8 


31 


36 


33 


2C9 


46 


22 


32 


2E1 


27 


16 


57 


4B1 


39 


59 


2 



Isotope Effects for ortho- and para-Methylanisole. Equal concentrations of 
p-MA and /^•^IA-a-^H3 were incubated with CYPs 1 A2, 2B1, 2B6, 2C9, 2E1 and 4B1 in 

15 order to measure intermolecular competitive KIEs. Identical experiments were carried 
out for the ortho isomers. Product formation by 2C8 was very low, most likely because 
a- and p-MA are poor substrates for 2C8; therefore, isotope effects could not be 
accurately measured for 2C8. For p-MA, a combined competitive KIE was calculated by 
dividing the isotope effect for /^-methoxybenzyl alcohol formation by the isotope effect 

20 for /?-cresol fomiation. The combined KIEs range from 5.9-6.8 with isoforms 1 A2, 2B1, 
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2B6, and 2E1, and the average is 6 ± 1 (Table 8). The values are 2,6 and 1 .4 for CYPs 
2C9 and 4B1, respectively. For o-MA, a combined KIE was calculated for each isoform 
by dividing the isotope effect for /^-methoxybenzyl alcohol formation by the product of 
the isotope effects for p-cresol formation and aromatic hydroxylation. The values range 
5 from 9 to 14, and the average is 1 1 ± 3 (Table 10). 

Intramolecular KIE were measured for o- and p-MA using o~ and 
/^■methylanisole-a-^Hl. respectively, with CYPs 1A2, 2B1, 2B6, 2C9, 2E1 and 4B1. The 
average KIE forp-MA was 4.4 ± 0.3, and the range of values was 3.69-6.08 (Table 8). 
For o-MA the average was 6.2 ± 0.6, and the range of values was 5.7-6.9 (Table 9). 
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Primary (P) and secondary (S) KIE were calculated for o- and p-MA by the 
method of Hanzlik et al. (X Am. Chem, Soc. 107, 7164-7, 1985) according to the 
equations 19-22. Here dl/dO is equal to the observed intramolecular KIE using o- or 
/7-methylanisole-a-^Hi as substrate (before statistical correction for the relative number 
of hydrogen to deuterium), and dO/d3 is the observed inteimolecular combined KIE. 

dl/dO = 2P/S (Eq. 19) 

dO/d3 = PS^ (Eq.20) 



S = 3p^ (Eq.21) 
V[(dl/d0)/2] 

10 P = [(dl/dO)/2]S (Eq. 22) 

Intramolecular KIE that are measured using dl substrates are assumed to be free 
from masking. Inteimolecular KIE measured with d3 substrates are free from masking 
when the value is greater than the intramolecular KIE value, therefore P and S were 

1 5 calculated only for isofonns that met this criteria. For p-MA, the range of P was 4.3-5,0, 
with an average was 5 ± 1 , and the range of S was 1 .08-1 .22, with an average of 1 . 15 ± 
0.22 for CYPs 1 A2, 2B1, 2B6, and 2E1 . For o-MA, the range of P was 6.7 -8.9, with an 
average was 8 ± 3, and the range of S was 1.15-1 ^9, with an average of 1 .2 ± 0.4 for 
isoforms 1A2, 2B1, 2B6, 2C9, 2E1, and 4B1. 

20 Total Product Formation with DO and D3 Substrates. Total product formation 

were independently measured for o- and /^-MA-a-^Ha and -a-^Ha in order to assess the 
amount of branching for both substrates. The data in Table 1 1 shows the amounts of 
MBA and p-cresol formed, expressed in imiol, from the metabolism of^MA-dO and d3 
for CYPs 1A2, 2B1, 2B6, 2C9, 2E1, and 4B1, and the products from /?-MA-dO substrate 

25 for 2C8. Metabolism of o-MA in similar experiments yields o-MBA, o-cresol, and a 

phenolic product identified as 4-methoxy-3-methylphenol (Table 12). For every enzyme 
and both of the isomeric substrates, the average amount of alcohol formed from d3 
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substrate is less than the amount formed from dO substrate, and the difference is 
statistically significant at the 0.05 significance level. For o- and p-MA, the average 
amount of cresol formed from 63 substrate is statistically greater than the amount formed 
from do substrate (a = 0.05), except when a-MA is metabolized by 2C9. The average 
5 amount of 4M3MPh formed from o-MA-dO is statistically different (a = 0,05) from the 
amount formed from o-MA-d3 with isofoims 1 A2, 2B1, and 2C9. No statistical 
difference was observed for CYPs 2B6, 2E1, or 4B1 (a = 0.05). 

Table 1 1 . Alcohol, Cresol and Total Product Formation from Metabolism of 

/>-Methylanisole and p-Methylanisole-a- H3 by Multiple CYP Isoforms 
C YP substrate nmol p-MBA (SD) nmol p-CR (SD) mnol total TP" (SD) " 



1A2 


DO 


23* 


(3) 


2* 


(4) 


25* 


(3) 




D3 


13* 


(2) 


8* 


(2) 


21* 


(3) 


2B1 


DO 


4.2* 


(3) 


2.0* 


(3) 


6.2* 


(5) 




D3 


0.9* 


(1) 


3.3* 


(3) 


4.2* 


(4) 


2B6 


DO 


7.0* 


(2) 


2.10* 


(2) 


9.1 


(2) 




D3 


2.39* 


(6) 


9* 


(2) 


11 


(2) 


2C8 


DO 


1.89 


(3) 


0.20 


(2) 


2.08 


(4) 




D3 


nd* 




nd 




nd 




2C9 


DO 


8* 


(1) 


0.42* 


(9) 


8* 


0) 




D3 


4.3* 


(4) 


1.16* 


(7) 


5.4* 


(4) 


2E1 


DO 


23.8* 


(9) 


5.6* 


(4) 


29.4 


(8) 




D3 


10.6* 


(7) 


18.2* 


(4) 


29 


(1) 


4B1 


DO 


63* 


(1) 


48* 


(!) 


111* 


(2) 




D3 


48* 


(1) 


52* 


(2) 


100* 


(1) 



1 0 Each result is the average of three, four or five individual determinations and the number 
in parenthesis indicates the standard deviation (SD) in the last significant digit in the 
value. " TP, nmole p-MBA + nmole p-CR. * Nd, not determined, due to low levels of 
product formation.. * Denotes statistical significance between the amount of product 
formed from metaboBsm of the DO substrate as compared to the D3 substrate (a = 0.05). 
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The total product (TP) fonnation is shown in Tables 1 1 and 12 for o- and p-MA 
as the sum of the individual amounts of products. The amount of TP formation from dO 
and d3 p-MA is statistically different at the 0.05 significance level for CYPs 1 A2, 2B1, 
2C9, and 4B1, and are not statistically different for 2B6 and 2E1 . Statistical comparison 
5 of TP from do and d3 o-MA (a = 0.05) showed differences for CYPs 2B1, 2C9, and 2E1, 
and no differences for 1 A2, 2B6, and 4B 1 . 

Regioselectivity of a-Chloromethyl-p-xy!ene. Metabolism of a-chloromethyl-^- 
xylene (a-ClMpX) by CYPs 1A2, 2B1. 2B6, 2C9, 2E1, and 4B1 yields a-chloromethyl- 
^benzyl alcohol and ^tolualdehyde (Figure 10). For CYP 2C8, the level of p- 

10 tolualdehyde that was formed was too low to be accurately measured. The 
regioselectivity results are listed in Table 13 as the ratio of metabolism at the 
a-chloromethyl group to metabolism at the benzylic position using a-ClMpX-dO as the 
substrate. The range of regioselectivity for 1A2, 2B1, and 2E1 is 1.8 - 5.4, which 
corresponds to a difference in a free energy of activation (AAG^) range of 0.3-1 kcal/mol, 

15 whereas for 4B1, 2C9, and 2B6 the regioselectivity range is from 0.035-0.76. 
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Isotope Effects for a-Chloromethyl-p-xylene, Combined intermolecular 
competitive KIEs were measured for a-ClMpX using equal concentrations of dO and 
d3 substrates. The combined KIE is the isotope effect for the formation of a- 
chloromethyl-/7-benzyl alcohol divided by the isotope effect for the formation of p- 
5 tolualdehyde. 

The combined KIEs range from 5.5-8.8 with isoforms 1A2, 2B1, 2B6, 2E1, 
and 4B1 (Table 9). The amounts of p-tolualdehyde formed by 2C8 and 2C9 were 
very low and could not be measured in the intermolecular competitive KIEs. 
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Intramolecular isotope effects were measured with a-chloromethyl-/>-xylene- 

a'-^Hi for seven isofoims (Table 13). The range is 6,0-8.7, and the average is 7 ± 2. 

The individual contributions from the primary and secondary isotope effects were 

calculated for 1A2, 2B1, and 2B6, since the intermolecular isotope effects for these 

5 isoforms were not masked. The average primary isotope effect is 7 ± 2, and the 

average secondary isotope effect is 1.0 ± 1. 

Total Product Formation with Do and D3 Substrates. The amoxmts of a- 

chloromethyl-j3-benzyl alcohol, /?-toluaIdehyde, and /7-toluic acid were measured and 

total product formation was calculated after the independent incubation of a-ClMpX- 

10 a'-*H3 and a-ClMpX-a'-^Hs with CYPs 1A2, 2B1, 236, 2C9, 2E1, and 4B1. p- 
Tolualdehyde was metabolized further by three isoforms 1A2, 2B1, and 2E1 to 
produce small amounts of p-Xo\mc acid. The amoimts of p-toluic acid formed and the 
percentages of the total amounts of product formed were 1 .3 ± 0.4 nmoles (3%) for 
1 A2; 0.7 ± 0.2 (3%) for 2B1 ; and 2.8 ± 0.2 (10%) for 2E1 . For all six isoforms, 

1 5 complete substitution of deuterium for hydrogen at the -CH3 group caused a decrease 
in benzylic hydroxy lation, when evaluated at the 0.05 level of significance, using the 
student's t-test (Table 14). 
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Table 14. Alcohol, Aldehyde, Toluic Acid and Total Product Formation from 
Metabolism of a-Chloromethyl-p-xylene and a-Chloromethyl-j^?- 



xylene-a'-^H3 by Multiple CYP Isoforms 



CYP 


substrate 


nmol D~Alc^ 


(SD) 


nmol p-Ald + 
TolA'iSV) 


nmolTP^(SD) 


1A2 


DO 


8,7* 


(3) 


41 


(4) 


50 


(4) 




D3 


1.3* 


(2) 


46 


0) 


47 


0) 


2B1 


DO 


7.9* 


(5) 


14* 


(2) 


22 


(2) 




D3 


1.3* 


(1) 


20.0* 


(4) 


21.3 


(5) 


2B6 


DO 


9,2* 


(5) 


7* 


(1) 


16* 


(1) 




D3 


1.82* 


(1) 


12* 


(1) 


14* 


(1) 


2C8 


DO 


0.11 


(1) 


nd^ 




nd 






D3 


nd 




nd 




nd 




2C9 


DO 


3.35* 


(3) 


0.4* 


(1) 


3.8* 


(1) 




D3 


0.10* 


(1) 


1.7* 


(3) 


1.8* 


(3) 


2E1 


DO 


4.5* 


(4) 


25 


(4) 


29 


(4) 




D3 


0.78* 


(3) 


30 


(1) 


31 


0) 


4B1 


DO 


62* 


(2) 


2.2* 


(2) 


65* 


(2) 




D3 


15.2* 


(5) 


1.87* 


(5) 


17.0* 


(4) 



Each result is the average of either three individual determinations, and the number in 
5 parenthesis indicates the standard deviation (SD) or propagated error (PE) in the last 
significant digit in the value. * a-Chioromethyl-/7-benzyl alcohol; **/?-Tolualdehyde; 
''/j-Toluic acid formed from the metabolism of j3-tolualdehyde; p-toluic acid 
formation was undetectable for CYPs 2B6, 2C9, and 4B 1 . TP, the total amount of 
products v^as calculated by the addition of the amoimts of individual products on this 
10 table. ^ Nd, not determined, due to low levels of product formation. *Denotes 

statistical significance between the amount of product formed from metabolism of the 
DO substrate as compared to the D3 substrate, using the student's t-test, at a 
significance level of 0.05. 
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For CYPs 2Bly 2B6, and 2C9, alternate product formation increased upon 
deuterium substitution whereas for 1 A2 and 2E1, the amount of alternate product 
formation with the -CD3 substrates is not significantly different from product 
formation with the -CH3 substrates (P = 0.05). There is a small, but statistically 
5 significant decrease in the amount of alternate product formation (a = 0.05) with 4B1 
when the substrate is deuterated. 

Discussion 

The structurzil/steric influence of an enzyme's architecture normally plays an 

1 0 important part in determining rates of product formation, regioselectivity and 

stereospecificity. When a large number of C YP isoforms are considered, it is easy to 
imagine that there is room for a great deal of variability in rates of product formation, 
regioselectivity and stereospecificity due to structural variability around the active 
site. This is true in many instances (Ortiz de Montellano, P. R. (1995) Cytochrome 

15 P450 structure, mechanism, and biochemistry, 2nd ed,. Plenum Press, New York). 
On the other hand, if regioselectivity is measured for a large number of isoforms, and 
the results are found to be the same for all en2ymes, it is possible that regioselectivity 
is controlled mostly by electronic factors, and the steric component is minor. If the 
steric component is minor, it is relatively certain that Kcq = 1 . 

20 Small substrates are likely to have fewer steric interactions with an enzyme, as 

compared to large substrates. It is likely that regioselectivity results for a single 
substrate would agree for a large number of CYP isoforms when steric factors are 
small or absent This statement is supported by isotope effect studies that argue for a 
conserved active oxygen for multiple CYPs. In addition, the relatively hydrophobic 

25 nature of the active site for most P450s reduces the possibility of strong electronic 
interactions between a small substrate and the enzyme. 

Evidence in the literature for the metabolism of the compounds shown in 
Figure 1 1 can be used to support the claim that steric factors are negligible or absent 
for small substrates with various CYPs. The metabolism of a group of small 

30 halogenated alkanes by CYP2E1 was described (both experimentally and 

computationally) by purely electronic factors (Yin et al. (1995) Froc. Natl. Acad Sci. 
USA 92, 1 1076-80). Iyer, et al. used isotope effects to show that o- and /^-xylene 
enjoyed rapid rotation in the active sites of CYPs from four sources, but that steric 



68 



BNSCOCID <WO 000e205A1_L> 



wo 00/08205 PCT/US99/1 771 3 

constraints prevented the larger 4,4'-dimethylbiphenyl molecule from moving freely 

in the active site (Tassaneeyakul et al. (1996) J. Pharmacol, Exp. Then 276^ 101-8). 

Finally, Jones et al. concluded that the large molecule benzo(12)pyrene and 

constraints in the active sites, possibly from the I-helix, of 15 different CYPs 

5 determine the stereospecific activation of the compound (Biochemistry 34, 6956-61, 

1995). The careful choice of small substrates in this invention allowed successful 

determination of which isoforms imposed active-site constraints on the substrates that 

were tested. 

The regioselectivity of o- and 7?-MA and a-chloromethyI-/7-xylene metabolism 
1 0 based on the series of isotope effect experiments in multiple CYP isoforms described 
above will now be discussed. 



P'Methylanisole 

The first step toward assessing whether or not rapid rotation of a substrate 

1 5 occurs in the active site of CYP is to compare intermolecular and intramolecular 

KIEs. For isoforms 1 A2, 2B1, 2B6, and 2EI, the intermolecular KIEs for /?-MA are 
greater than the primary KIEs, which were calculated from the intramolecular KIE 
values. The fact that the intermolecular KIEs are greater than the corresponding P 
values is not unexpected since the intermolecular isotope effects are a combination of 

20 a primary and two secondary isotope effects. The rule of the geometric mean states 
that isotope effects at a given reaction center are independent and multiplicative, and 
the data is consistent with this rule. Therefore the inter and intramolecular KIEs are 
equal for these isoforms. 

A comparison of the amounts of total product formation for Do and D3 

25 substrates was used to assess the effect of substrate deuteration on branched 

pathways. Upon deuteration of the substrate at the benzylic position, the amount of 
/7-MBA is reduced, and the amotmt of p-cresol is increased, and the changes are 
statistically significant at the 0.05 significance level, for all isoforms tested (Table 
1 1). Total product formation for Do and D3 substrates with 1A2 and 2B1 are 

30 statistically different (a = 0.05). The small differences in total product formation can 
be best explained by the fact that branching to water occurs. If rapid rotation occurs 
between the ES complexes, formation of a second alternate product only serves to 
unmask reactivity differences. 
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A number of results uphold the concept of rapid "rotation" of some substrates 
in the active site of CYP and no change in overall product formation upon 
deuteration. Lindsay Smith et al. performed experiments with anisole (Journal of the 
Chemical Society Perkins Transactions II, 621-8, 1983). Although they did not 
5 determine the intramolecular isotope effect they found that deuteration of the methyl 
group resulted in significant branching to para hydroxylation of anisole. This 
resulted in an intermolecular isotope effect of 7. This indicates that the ratios 
observed by these investigators might reflect the intrinsic difference in reactivity of 
the aromatic and demethylation pathways. Lindsay Smith et al. observed no 
10 significant change in overall turnover for the and Do anisole. Thus, all the 

branching gave the alternate product, 4-methoxyphenol, and water formation was not 
responsible for the large isotope effect. The relative difference in energy for p- 
hydroxylation of anisole and demethylation can be calculated to firom the product 
ratios to be 0.6 kcal/moL 
15 Another example of efficient branching to an alternate position is the 

metabolism of 7-ethoxycoumarin by CYP1A2. Deuteration of the ethoxy group in the 
methylene position gave rise to a decrease in the deethylation pathway, an increase in 
aromatic hydroxylation, and no change in the overall product formation or NADPH 
consumption. The energy difference between the two pathways was determined to be 
20 about 1 .8 kcal/mol, with O-deethylation as the favored path. The above data for 1 A2, 
2B1, 2B6, and 2E1 is consistent with these data found in the literature in that rapid 
rotation of substrate occurs in the active sites, and intermolecular KIEs are unmasked. 

The data were evaluated in terms of criteria 1 and 2. Since the intermolecular 
KIE > primary KJE, and since the total product formation does not not significantly 
25 changed upon deuteration of the substrate, it is concluded that rotation of the substrate 
in the active sites of CYPs 1 A2, 2B1, 2B6, and 2E1 is rapid. The regioselectivity 
values for isoforms 1A2, 2B1, 2B6, and 2E1 are similar, which indicates that Keq = 1. 
Together this data indicates that the observed regioselectivity is not influenced by 
steric constraints fi-om the enzyme active site, and that benzylic hydroxylation is 
30 energetically more favorable than O-dealkylation by 0.5-1 .8 kcal/moL 

The intermolecular and primary KIEs for benzylic hydroxylation of p-MA by 
CYP 2C9 are not equal. The most likely explanation for this observation is that the 
amount of branching of the metabolism to the alternate product is very low, which is 
indicated by the extremely high observed regioselectivity of 20:1 for BzH.ODAlk. 
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The regioselectivity suggests a high k^lh ratio (see Eq. 18) which can mask the KIEs 

and regioselectivity differences. 

The inter- and intramolecular KIEs for benzylic hydroxylation of /?-MA by 

CYP 4B 1 are not equal. For 4B1, it is not likely that the amount of branching is 

5 insufficient to unmask the isotope effect since the regioselectivity, expressed as 

BzH:ODAIk is very low, 1 .3; in fact 4B 1 shows the highest amoimt of ODAlk relative 

to BzH for the series of isoforms that were tested. The reason for the nonequivalent 

isotope effects can be explained either by slow rotation of the substrate in the active 

site or steric effects from the enzyme that causes Keq ^ 1 (with or without rapid 

10 rotation). 



o-Methylanisole 

The intermolecular and intramolecular KIE for o-MA with CYPs 1 A2, 2B 1, 
2B6, 2C9, 2E1, and 4B1 are equal, which support rapid rotation of the substrate in the 

15 active sites of the enzymes. The regioselectivity, expressed as BzHrODAlk for 1A2, 
2B1, 2C9, and 2E1 ranges from 1.2-2.9, while the regioselectivity for 2B6 is 0.52 and 
for 4B1 is 0.7, which are clearly outliers. When expressed as B2H:Arom, the 
regioselectivity for 1A2, 2B1, and 2E1 are in agreement (the range is 0.3-0.6), the 
value for 2B6 is 2 and for 4B1 is 22. The value for 2C9 has a large standard error and 

20 cannot be accurately assessed in relation to the other isoforms. The results suggest 
that Kcq = 1 for o-MA with IA2, 2B1, and 2E1 but not with 2B6 or 4B1 . Thus steric 
constraints on o-MA due to the enzyme active site are apparently absent in 1 A2, 2B1, 
and 2E1, and the rank order of reactivity is aromatic hydroxylation > benzylic 
hydroxylation > O-dealkylation. This data agrees with the data for p-MA and CYPs 

25 1 A2, 2B1, 2B6, and 2E1 in that benzylic hydroxylation is energetically more 
favorable than O-dealkylation. 

When the ratio of branching increases relative to metabolism at the benzylic 
position with the ortho isomer, the inter- and intramolecular KIE for 2C9 are 
unmasked, and the regioselectivity results agree with the majority of the data for the 

30 other isoforms in that BzH is preferred over ODAlk. Deuteration of the substrate 
does not cause a significant increase in alternate product c?-cresol, and the average 
amount of 4M3MPh is 2 for dO substrate, and 1.8 for d3 substrate, a small or 
negligible change. Total product formation for o-MA-d3 is approximately half of the 
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value for the dO substrate. It is possible that branching to water is responsible for 

unmasking the intermolecular KIE for o-MA and 2C9, 

Inspection of the data for o-MA with 4B 1 shows interesting results. The 

intermolecular isotope effect is not masked and total product formation for dO- and 

5 d3 -substrates is equal. However, the regioselectivity is the opposite of the values 

observed for the majority of isoforms and substrates. It is probable that rapid rotation 

of the two functional groups in the active site occurs since the benzylic and methoxy 

groups are in close proximity to each other. Therefore, it is not surprising that the 

intermolecular isotope effect is unmasked for 4B 1 . The fact that the amount of 

1 0 phenol product from o-MA is extremely low together with the inverse regioselectivity 

results, expressed as BzHiODAlk, suggest that Keq ^ 1 or, in terms of the Curtin- 

Hammett principle, the transition-state for benzylic hydroxylation is sterically 

hindered. 

It is proposed as part of the invention that the active site of 4B1 may have a 

15 narrow channel over the heme, which does not permit fast rotation of the substrate 
para-MAf but which allows the p-MA to translate away from the heme, flip over, and 
translate back towards the heme. The proposed narrow channel can accommodate 
either of the two functional groups, but accommodates the methoxy group better. 
This is indicated by (i) the fact that regioselectivity for p-hAA shows the lowest value 

20 when compared to the majority of isoforms tested, and (ii) metabolism of the methoxy 
group is preferred over metabolism of the benzylic group for o-MA, which is opposite 
to the results for the majority of isoforms. Almost no ring hydroxylation is observed 
when o-MA is metabolized by 4B1 , in contrast to other isoforms. This also supports 
the idea of a preferred binding mode. 

25 It has been suggested that the presence of an amino, substituted amino group, 

or an oxygen atom may confer substrate specificity for 4B1 (80). The presence of a 
positively charged amino acid in the active site of 4B1 , near the heme, may confer an 
electrostatic component to the binding and regioselectivity. A positive charge could 
interact with either N or O in the substrates, and bias the position of the substrate in 

30 the active site. 

The isotope effect results for o-MA with 2B6 suggest that the substrates 
undergo rapid rotation in the active site. However, the regioselectivity values for 
2B6, when expressed as BzH:ODAlk and BzH:Arom, are opposite to the 
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regioselectivity values for 1 A2, 2B1, and 2E1 . Data in the literature agrees with the 

results disclosed herein. Metabolism of toluene by expressed 2B6 gave a high ratio of 

BzH:ODAlk when compared to other isoforms. It is proposed as part of the invention 

that the active site of 23 6 is large enough to allow for the rapid rotation ofpara-MA 

5 and that 2B6 has a narrow channel in the active site near the heme. The narrow 

channel is responsible for constraining the ortho isomer of MA in one orientation. In 

conclusion, the present data suggests that K«q ^ 1 for the 0-MA-2B6 complex. 

An additional observation of the present regioselectivity data for o-MA is that 

the ratios of ring hydroxy lation para to the methoxy group to ODAlk, according to 

10 isoform, are 4.5(1A2), 3.5(2B1), 1.5(2C9), and 3.5(2E1). These data are similar to 
the ratios of p-hydroxyanisole to phenol (average ratio = 2.9) that were observed 
when anisole (Ph-OCHa) was metabolized by four microsomal systems. The present 
data and Lindsay Smith's data are in contrast to the ratios of ODAlk :Arom measured 
herein for 2B6 and 4B1, which are 0.3 and 0.03, respectively. This lends additional 

1 5 support to the hypothesis that steric factors contribute significantly to the observed 
regioselectivity and KIEs for o-MA with 2B6 and 4B 1 . 

The degree of symmetry in the transition-state can be assessed using the P 
values calculated from the intramolecular KIEs for metabolism at the ben2ylic 
position. The average P value for ;7-MA with CYPs 1A2, 2B1, 2B6, and 2E1 is 5 ± 1, 

20 and the average P for <?-MA is 8 ± 3. The average P value for octane was 7.6, which 
was postulated to have a symmetrical transition-state (Jones et al. (1987) 7. Am. 
Chem. Soc. 109^ 2171-3). The average P value for p-MA of 5 implies an asymmetric 
transition-state and an exothermic reaction relative to H-atom abstraction from octane. 
The average P value for o-MA suggests a symmetrical transition state. It is suggested 

25 as part of the invention that the /^-methoxy group provides resonance stabilization in 
the transition-state and accounts for the lower P value for p-MA as compared to the P 
value for octane (see Figure 12). 

The reason for the differences in the average intramolecular KIEs for o- and 
p-MA is not known, but is most likely explained by the fact that the two isomers are 

30 electronically different. According to the Melander and Westheimer (Melander, L. 
(1 960) Isotope Effects on Reaction Rates ^ Ronald Press, New York), the magnitude of 
an intrinsic KIE is governed by the symmetry of the transition-state. Therefore, the 
present data indicates that the transition-state for the carbon-hydrogen bond breakage 
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of/?-MA is earlier than the transition-state for o-MA, Based on the Hammond 

postulate, the reaction for p-MA is more exothermic, and the product is more stable, 

when compared to a-MA. This comparison is based on the assumption that the 

transition states for both isomers are linear, which may not be true since steric 

5 hindrance in the ortho compound could induce a non-linear transition state. But the 

effect of a non-linear transition-state would be expected to lower an isotope effect. If 

the para compound is used as a reference (and was linear) the ortho isomer would be 

expect to have a lower intrinsic KIE as compared to the para isomer. It is possible 

that resonance effects may be responsible for differences in stabilization of the ortho 

1 0 and para reactants and products, and thus the different KIBs. 

a-Chloromethyl-p-xylene 

The results in Table 1 3 show that the intermolecular KIE = primary KIE. 
Rapid rotation of the substrate in the active site is inferred for CYPs 1 A2, 2B1, and 

15 2B6. Among these three isoforms, 1A2 and 2B1 show similar regioselectivity results. 
Total product formation is not decreased upon deuteration of substrate for 1A2 and 
2B 1 ; therefore water formation is not responsible for unmasking the KIE. For 1 A2 
and 2B6, it is concluded that K«, = 1, and that metabolism at the a-chloro-methyl 
group is energetically more favorable than at the benzylic position by 036-0.93 

20 kcal/mol. 

The preference for metabolism at the a-chloro-methyl group, as compared to 
the benzylic position, is the predicted preference based solely on a comparison of the 
stability of the resultant radical products of hydrogen atom abstraction. The product 
formed from hydrogen atom abstraction at the carbon atom adjacent to the CI is 
25 resonance stabilized by CI atom. Weaker, inductive effects from the para-positioned 
benzylic group stabilize the product formed from hydrogen atom abstraction at the - 
CH3 group. 

While the isotope effect results for 2B6 indicate rapid rotation of the substrate 
in the active site, the regioselectivity data suggest that steric constraints in the active 
30 site are significant. This is based on comparison of regioselectivity results to the 

results for 1 A2 and 2B1. The results also correlate with the results for interaction of 
o-MA with 2B6, which suggested binding-site restrictions. 

The regioselectivity data for 2C9 cannot be assessed for rapid rotation, but an 
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interesting point can be made based upon a comparison of the results to the observed 

regioselectivity for p-MA. Metabolism at the -CH3 group of />-methyiaiusole and 

a-chioromethyl-p-xylene by 2C9 is highly favored over metabolism at the positions 

para to this methyl group, by 20:1 and 8:1, respectively, when compared to most of 

5 the other isoforms that were tested. It is possible that the proposed region of positive 

charge in the active site of 2C9 (Jones et al. (1996) DrugMetab. Dispos. 24^ 1-6) 

interacts with the lone pair(s) of electrons from the oxygen in the methoxy group and 

the chlorine in the a-chloromethyl group. This would act to keep the substrates 

contained, and subsequently influence the regioselectivity results by this electrostatic 

10 interaction. 

The intermolecular KIE is not equal to the intramolecular KIE for 2E1, which 
indicates that rotation of the substrate is not fast in the active site. Incubation of the 
ortho isomer of this compound would potentially permit rapid exchange of the 
substituents in the zictive site, and would allow the assessment of regioselectivity 

15 results without steric or other constraints. 

The isotope effect results for 4B1 signify that rotation of the substrate is not 
rapid. The regioselectivity data shows a drastic preference for benzylic hydroxy lation 
to metabolism at the a-chloromethyl group, the preference is -30: 1 . Deuteration of 
the substrate shows a statistically significant, but very small, decrease in alternate 

20 product formation, from 2.2 nmol to 1 .87 nmol. The regioselectivity and 

noncompetitive product data suggest a large structural constraint imposed by the 
enzyme. Another possibility for potential 4B1 -substrate interaction is the presence of 
a positively charged region in the active site. This idea was used to explain the 
regioselectivity data for o-MA, and 4B1 substrates preferences (Blaise Smith et al. 

25 (1995) Biochem. Pharmacol 50, 1567-75). A list of known 4B1 substrates contains 
compounds that usually have an amino group, a substituted amino group, or an 
oxygen atom. 

The accurate determination of whether or not Kcq = 1 (criterion 2) is the most 
30 challenging element of the use of the disclosed methodology for separating electronic 
and steric factors. The present invention's method of comparing regioselectivity 
values for multiple isoforms is an indirect method. While the invention concludes 
that the concept is valid, there are some potential shortfalls. The technique requires 
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the use of a large number of isoforms for metabolism studies. If compounds are not 
good substrates for a large enough number of isoforms, the sample size for assessing 
criterion 2 will be small. If regioselectivity veilues vary widely for a small sample 
size, it may not possible determine which isoform(s), if any, impose steric constraints 
5 on the substrate. 

Among the present data for a-chloromethyl-/7-xylene, only three isoforms 
(1 A2, 2B1, and 2B6) show rapid rotation of the substrate in the active site, as judged 
by the comparison of intermolecular and intramolecular isotope effect values. 
Regioselectivity values for two of the three isoforms were in agreement (1 A2 and 

10 2B 1 ), and the other isoform (2B6) showed the opposite regioselectivity value. With a 
sample size of only three, it is difficult to prove that the two isoforms whose values 
agree are not outliers, and that the distinct value does not represent a larger set of 
consensus data. However, information was gained by the comparison of groups of 
data from multiple substrates. For instance, it was concluded that isoforms 2B6 and 

15 4B1 imposed structural constraints on the substrate o-methylanisole, and the data for 
a-chloromethyl-/?-xylene was evaluated with this information in mind. For both o- 
methylanisole and a-chloromethyl-/^xylene, the regioselectivity data for isoforms 
2B6 and 4B1 agree, and in both cases the regioselectivity values indicate that 
metabolism occurs more often at the less energetically favorable site (judging from 

20 the predicted stability of ±e proposed carbon based radical product). The 

consistency of the data lends support to the invention. Regioselectivity values and 
isotope effect experiments Avith additional isoforms as well as additional substrates 
can be obtained in order to strengthen the initial hypotheses. 

From these results, as well as the results of Lindsay Smith et al. (Journal of 

25 the Chemical Society Perkins Transactions II, 621-8, 1983) and Harada et al, (J. Biol. 
Chem. 259, 3005-10, 1984), a rank order of reactivity of functional groups that are 
metabolized with freedom from structural constraints from the apoprotein can be 
established. The following order of reactivity: a-chloromethyl hydroxylation > 
benzylic hydroxylation > 0-deethylation > aromatic oxidation > O-dealkylation is 

30 predicted. While Lindsay Smith et al. and Harada et al. did not perform a rigorous set 
of isotope effect experiments to check for rapid rotation of substrate, stoichiometric 
results from both data sets indicated that rapid rotation was likely. 

In simimary, a methodology for determining whether the electronic component 
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of CYP -mediated reactions are free frota the structural component has been 

developed. The relatively small size of these substrates is partially responsible for the 

success of this model. Electronic models correctly predict the primary site of 

metabolism for about 80% of substrates. For CYP enzymes that do not permit rapid 

5 rotation of substrates, an electronic component still exist in combination with the 

steric component. As such, electronic models are still required in predictive models 

for most CYP enzymes. 

The regioselectivity results from these experimental studies can be compared 

to, and used for the validation of existing computational or chemical methods for 

1 0 predicting CYP reactions. In addition, the experimental data can be used to help 
refine computational models, and as guidelines for the development of new models. 
Ultimately, methods for predicting the electronic component to reactivity can be 
combined with methods for predicting the influence of protein structure on 
metabolism. Ideally, this will allow for the construction of complete, more accurate 

1 5 models, that is, models that will incorporate the influence of both electronics and 
steric factors, without imposing limits on substrate size. 

In this example, chemicals and reagents were obtained as described in 
Example 2. 

Synthesis of Alcohols and Toluenes. Methyl-4-methoxybenzoate and methyl- 
20 2-methoxybenzoate were used as starting materials for synthesis of the corresponding 
substituted benzyl-a- H2 alcohols, and substituted toluenes-a, a, a- H3 and -a- Hj 
were synthesized from the corresponding D2 and Do alcohols, respectively, and 
purified according to the procedures described in Examples 2 and 3 

Synthesis of 4'Methoxy-3-methyl-phenol. The synthesis for the phenol was a 
25 two step synthesis via the Vilsmer-Haack condensation to obtain the aldehyde and 
then by the Baeyer-Villiger oxidation to give the phenol. Phosphorus oxychloride 
(5.1 mL, 0.057 mol) was added the dropwise to a solution of 2-methylanisole (5.33 
mL, 0.043 mol) in dimethylformamide (6 g) under nitrogen. After addition, the 
mixture was heated and refluxed for 4 h, cooled then added to 1 00 mL of water. An 
30 excess of 1 0% NaOH was added and the solution was extracted with 4 x 1 00 mL of 
ether. The ether layers were washed with water and then brine solution, dried over 
magnesium sulfate and evaporated under reduced pressure. The dark oil that 
remained (5.0 g, 77%) was purified by column chromatography (silica gel: 90% 
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hexane, 10% ethyl acetate). The remaining oil (1.78 g, 28%) was dissolved in 
dichloromethane (100 mL) and then 3.03 g, (0.018 mol) 3-chloroperbenzoic acid was 
added. The mixture was refluxed under nitrogen for -5 h, cooled, and then 
evaporated under reduced pressure. The residue was dissolved in ethyl acetate and 
5 washed with aqueous saturated sodium bicarbonate. The organic layer was dried over 
magnesium sulfate and evaporated under reduced pressure. The dark oil that 
remained was dissolved in 7 mL of methanol and 5.5 mL of 10% aqueous potassium 
hydroxide solution was added. The solution was stirred under nitrogen for ~2 h, 
water was added, and the solution was extracted with ether. The ether layers were 
1 0 combined and washed with saturated sodium bicarbonate solution and then with 
water. The ether was dried over magnesium sulfate and evaporated under reduced 
pressure. The oil that remained (1 ,21 g, 74%) was recrystallized twice with hexane to 
give long white needles (0.34 g, <1%). The purity was assessed by NMR and GC- 
MS. 

15 Synthesis of a-Chloromethyl-p-bemyl-CL- H2 Alcohol To a solution of 4- 

(chloromethyl)benzoic acid (2.5 g, 0.01 5mol) in ether was added an excess of 
diazomethane in ether. After the excess of diazomethane had evaporated from the 
reaction, the ether was evaporated under reduced pressure leaving solid (2.68 g, 99% 
crude yield) which was then dissolved in ether and added dropwise to a solution of 

20 lithium al\iminum deuteride (0.73 g, 0.017 mol) in ether cooled to 0 °C. A wliite 
semi-solid was isolated (2.3 1 g, 99%) after work-up as per the para-substituted 
benzyl-a-^Hi alcohols in Example 2. This compound was used without further 
purification. 

Synthesis of a-Chloromethyl-p-benzyl AlcohoL The reaction was performed 
25 using the same method as the method described above for the synthesis of a- 

chloromethyI-7>-benzyl-a-^H2 alcohol, except lithium aluminum hydride was replaced 
with lithium aluminum deuteride (82% yield). 

Synthesis of a-Chloromethyl-p-xylene-dS and a-Chloromethyl-p-xylene-dL 
These were prepared in the same manner as the methylanisoles. The oils that 
30 remained were first purified by silical gel column chromatography (1 00% hexane) 
and then an alumina column (100% hexane). 

Preparation ofHepG2 Expressed CYP, Incubation Conditions, Controls, 
Product Isolation, and Data Analysis. The methods and details of protein expression 
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and preparation for isoforms 2B6, 2B8j and 4B1, incubation conditions, sample work- 
up, and data analysis were the same as described above in Example 2 except for the 
following details. Competitive isotope effect experiments were performed in 
triplicate with mixtures of Do and D3 substrates, each at a final concentration of 0.125 
5 \iM. The independent measurements of product formation from Do and D3 substrates 
were performed at final concentrations of 0.25 p,M. Product formation from Do 
substrates was measured with 10 nmol of the corresponding substituted benzyl 
alcohol dideuterated at the benzylic carbon as the internal standard; for D3 substrates 
the 10 nmol of the corresponding nondeuterated benzyl alcohol internal was used as 

10 the standard. o-Cresol,/7-cresol, and 4-methoxy-3-methylphenol were quantified with 
20 nmol of /?-bromobenzaIdehyde as the internal standard, against a standard curve. 
A small amount of background o- or /?-cresol was quantified and subtracted fi-om the 
observed amount of the corresponding cresol product 

Product Analysis by GC-MS, An HP-1 column was used for product 

1 5 quantification with the same GC-MS run conditions as those stated in the 

Experimental section of Chapter II except the oven conditions were 50 °C for 0.5 min, 
1 0 °C/min to 145 °C and 20 °C /min to 250 °C, and 2 min at 250 °C. An HPWax 
column was used for the identification of 4-methoxy-3-methylphenoI. Refer to 
Example 2 for the remaining details. 

20 Identification of o-Methylanisole Metabolites. Three product peaks were 

observed after GC-MS analysis of sgimples from incubations of o-methylanisole and 
CYP. Two of the peaks were imambiguously identified as 4-methoxybenzyl alcohol 
and o-cresol by comparison to authentic standards. A third peak for all isoforms 
exhibited mass spectral fragment patterns and GC retention times on both HP*1 and 

25 HPWax columns that were identical to the authentic compound 4-methoxy-3- 

methylphenol. Since a mass spectral fragment pattern is normally unique for any 
given chemical, the third peak has been labeled 4-methoxy-3-methylphenol. 

However, authentic compounds corresponding to the three other possible 
phenolic regioisomers were not available for comparison to the third peak. If the 

30 mass spectral fragment patterns, retention times, and area coimts relative to internal 
standard are the same for all four isomers, then the data for 4-methoxy-3- 
methylphenol represents a sum of the ring hydroxylation products. All o-MA 
incubation samples were analyzed on an HPWax column, which is more polar than an 
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HP- 1 column. The HPWax column was used to separate peaks for other isomeric 

compounds that normally overlap on an HP-1 column. Analysis of all o-MA 

incubation samples on the HPWax column did not result in any degree of peak 

separation. Therefore, it is reasonably certainty that the third product is 4-methoxy-3- 

5 methylphenol. 

Statistical Analysis, All statistical analysis was performed with the student's t- 
test assuming equal variances xising the Data Analysis Toolpak in Microsoft Excel 97. 

10 Example 5 

Regioselectivity as a measure of relative rate constants for hydrogen atom abstraction: 

As explained above, experimental metabolite ratios may be used to 
parameterize the theoretical model for CYP reactions given rapid rotation of the 
substrate, indicating the absence of enzyme imposed steric constraints. An additional 
15 means of confirming such rapid rotation in CYP enzymes is by examination of the 
relative rates of benzylic hydroxylation and 0-dealkylation of 4-methylanisole, The 
observed ratio of products from the two reactions permits confirmation of whether the 
reactions can be simply modeled with electronic factors. 

Table 1 5 demonstrates that for reactions on 4-methylanisole by various CYP 
20 enzymes, rapid rotation rates are present. 



Table 15. The regioselectivity, and inter- and Intramolecular kinetic isotope 

effects for 4-methylanisole 


Enzyme 


Benzylic/ 
0-Deaif^Iation 


(kcal/mole) 


Intramolecular 
KIE 


Intermolecular 
KIE 


CYP2E1 


5.41 


1.06 


4.02 


6.53 


CYP3A4 


12.83 


1.61 


5.27 


2.76 


CYP2B6 


3,93 


0.86 


4.58 


7.60 


CYP2C8 


8.26 


1.32 


3.60 


3.60 


CYP2C9 


19.62 


187 


3.92 


3.52 


CYP1A1 


5.91 


1.12 


4.23 


2.17 


CYP1A2 


17.98 


1.82 


3.79 


5.94 


Mean 


10.60 (5.85) 


1,38 (0,39) 


4.20 


4.59 



In addition to the normal isotope effects observed for the benzylic 
hydroxylation reactions, inverse isotope effects are observed for demethylation. This 
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indicates that the rate of interchange of the methyl and methoxy groups is fast for 
most en2:ymes tested. All enzymes favor benzylic oxidation, with an average different 
in AG^ of 1 .3 kcal/mole. This standard deviation is within the predictive error of the 
models. Thus, this data further supports the use of experimental metabolite ratios to 
5 provide relative rate constants for oxidation of different functional groups. 

To parameterize hydrogen atom abstraction, the results of a recent study 
proved exciting, not only in terms of mechanistic insight but because it presented with 
a potential apoprotein independent model for CYP mediated reactions. AVhile 

1 0 studying the mechanism of N-dealkylation reactions, a number of experiments were 
conducted with the well characterized hydrogen atom abstractor t-butoxy radical. 
These studies extended previous work by others which uses a laser to flash photolyze 
t-butylperoxide to t-butoxy radical. The classic method uses diphenylmethanol to 
follow the rate of hydrogen atom abstraction of a competing reaction. 

15 Of particular interest was comparing the isotope effect profiles for t-butoxy 

radical hydrogen atom abstraction, the deprotonation of a radical cation after outer- 
sphere electron transfer and CYP mediated reactions. It was found that not only were 
the isotope effect profiles similar, they were essentially identiczil for the t-butoxy and 
CYP mediated reactions with substituted dimethylanilines. This led to the hypothesis 

20 that t-butoxy radical is a good model for CYP mediated reactions in general. 

Using BLYP/6-31G* density function methods to obtain ground state 
descriptors the linear relationship (see Equation 24 and Figure 14) between 14 
experimentally measured rate constants (In rate) for hydrogen atom abstraction by 
t-butoxy radical and the BLYP/6-3 IG* heat of reaction (Hrcac) and radical 

25 intermediate HOMO energy in kcal/mol (Ehomo) was obtained. Hie standard error of 
the estimate for this correlation is aroimd 0.9 kcal/mole (see Figure 14). 

In rate = 12.43 + 0.22 (AHr) + 1 .84 (Ehomo) Eq. 24 

30 These results support the use of chemical systems to experimentally 

parameterize computational models. 
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Example 6 

Regioselectivity as a measure of relative rate constants for aromatic oxidation: 

Experimental metabolite ratios may also be used to parameterize the aromatic 
oxidation reaction in a theoretical model for CYP reactions. As with hydrogen atom 
5 abstraction discussed in the previous Examples, experimental ratios of aromatic 

oxidation to an alternate position may be determined using various CYP enzymes. An 
example of substrates that can be used to parameterize aromatic hydroxylation 
reactions are given in Figure 13. These substrates were chosen based on their small 
size, the presence of an alternate position for metabolism, and their hydrophobicity. 

10 In order for the substrates to be useful, the average of the rotation rates in the 

active sites should be rapid. For five of the substrates, 7-ethoxycoumarin, toluene, 
o-xylene, p-xylene, and anisole, isotope effect studies have shown that rapid exchange 
of the relevant positions does occur. The remainder of the substrates cannot form 
hydrogen bonds. Therefore, hydrophobic interactions are most likely to be involved 

15 in binding of these substrates, and rotation is likely to be rapid. Several human 

expressed CYP enzymes were used in these studies, including CYPl Al, 1 A2, 2A6, 
2B6, 2C8, 2C9, 2D6, 2E1, 3A4, and 3A5. The use of multiple enzymes insures that a 
reasonable average of any orientation preferences within the active site (see Figure 1). 
The analysis and construction of a model for both hydrogen abstraction and 

20 aromatic oxidation will require a combination of data from both reactions. This is 
accomplished using Equations 5 and 6 above. 

All references cited herein are hereby incorporated by reference in their 
entireties, whether previously specifically incorporated or not. 
25 Having now fully described this invention, it will be appreciated by those 

skilled in the art that the same can be performed within a wide range of equivalent 
parameters, concentrations, and conditions without departing from the spirit and scope 
of the invention and without undue experimentation. 

While this invention has been described in connection with specific 
30 embodiments thereof, it will be understood that it is capable of further modifications. 
This application is intended to cover any variations, uses, or adaptations of the 
invention following, in general, the principles of the invention and including such 
departures from the present disclosure as come within known or customary practice 
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within the art to which the invention pertains and as may be applied to the essential 
features hereinbefore set forth as follows in the scope of the appended claims. 
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1. A method of detennining the electronic tendency for metabolism of a 
compound by an enzyme with broad substrate specificity comprising the steps of 
5 a) determining, without reference to said enzyme, the regiospecific reactivity 

of each functional group in said compound towards metabolic transformation by any 
reaction(s) mediated by said enzyme; 

b) deriving one or more equations relating the activation energy for each 
enzyme mediated reaction to quantum chemical descriptors; 
10 c) determining, from said one or more equations, the metabolism of said 

compound by predicting the activation energies for each reaction; 

wherein said determined electronic tendency provides the electronic 
component for said reactions. 



15 

2. The method of claim 1 wherein step b) includes relating the activation 
energy for each enzyme mediated reaction to the heat of reaction with said compound. 

3. The method of claim 1 wherein step b) includes basing said one or more 
20 equations on semiempirical methods. 

4. The method of any one of claims 1-3 wherein step b) fiirther includes 
refining said one or more equations by parameterizing with experimental results to 
correct errors in activation energies. 

25 

5. The method of claim 4 wherein said experimental parameterization 
includes identifying the relative rates of substrate rotation in the active site of the 
enzyme by determining the isotope effect profiles for said reaction(s). 

30 6. The method of any one of claims 1-3 wherein step b) further includes 

refming said one or more equations by calculating high level quantum chemical 
parameters to correct errors in activation energies. 
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7. The method of claim 6 wherein step b) further includes application of ab 
initio methods to correct the quantum chemical activation energy (AMI) and/or the 
spin contamination in said one or more equations. 

5 8. The method of any one of claims 1-3 wherein step c) comprises identifying 

whether metabolism of the compoimd results in one or more of the following: 

(1) the production of an excess concentration of a toxic metabolite; 

(2) an excessive rate of metabolism of a pharmaceutically usefid substance; 

(3) an inadequate rate of metabolism of a pharmaceutically useful substance; 

10 and 

(4) an excessive rate of metabolism of a non-pharmaceutically useful 
substance into a toxic metabolite. 

9. The method of any one of claims 1 -3 wherein said any enzyme mediated 
1 5 reaction(s) includes hydrogen atom abstraction and/or oxygen addition. 

10. The method of claim 9 wherein said reaction is hydroxylation and/or 
aromatic oxidation. 

20 11. The method of claim 1 0 wherein said enzyme is a monooxygenase. 

12. The method of claim 1 1 wherein said monooxygenase is cytochrome P- 

450. 

25 13. The method of claim 12 wherein said cytochrome P-450 is selected from 

the group consisting of human enzymes CYP2E1, CYP3A4, CYP2B6, CYP2C8, 
CYP2C9, CYPlAl, CYP1A2, CYP2C19, CYP2D6, CYPIBI, and CYP2A6. 

14. A method according to claim 1 wherein said reactions are sterically driven 
30 with an electronic component. 
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15. The method according to claim 1 wherein said functional group is selected 

from the group consisting of C-H; C-C; C=C; C=C; C=0; C-N; C=N; -S-; -N-; -N=; - 

CHO;-OH;and-COH. 
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Figure 1. A proposed catalj^c cycle for cytochrome P450 (CYP). 




Phenylacetonitrile 




7 

Benzo[a]pyrene 



1 0 Figure 2. Cytochrome P450 (CYP) substrates which can be metabolized to 

toxic compounds. 
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9 I ! I I . 1 . ! ^ 1 ^ i ' ' ' 1 ' ' 

9 10 11 12 13 14 15 16 17 
Calculated Activation Energy 

Figure 3 depicts results of hydrogen abstraction reactions with the methoxy radical. 
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Methoxy Addition to Aromatic 
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5 Figure 4 shows correction of AMI heats of reaction to match DFT values 
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Figure 5 shows experimental AAG versus predicted AAH for ratios of hydrogen 
abstraction and aromatic oxidation reactions 
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Figure 6. Isotope effect profiles for six para-substituted toluenes with expressed 
5 CYPs. 
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OCH3 

w '''' 

p-Methylamsole 

o-Methylanisole 





H2CIC <f CH3 



a-Chloroniethyl-/?-xylene 

Figure 7. Cytochrome P450 (CYP) substrates used for the measurement of 
regioselectivity. 
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OCH3 




P450 + 




CH3 



CH3 



OCH3 




CH20H 

10 Figure 8. The formation of /7-cresol and p-methoxy benzyl alcohol from /?- 

methyianisole by CYP. 
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p-methylanisole 



HOH2C 




jEJ-methoxybenzyl alcohol p-cresol 

Figure 9. Two metabolites from the CYP mediated metabolism of 
5 /7-methylanisole. 




a-chloromethyl-/?-xylene 




a-chloromethyl-/>-benzyl alcohol tolualdehyde 

Figure 10. Two metabolites from the CYP mediated metabolism of a- 
chloromethy l-/>-xylene . 
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o-Xylene />-Xylene 




4,4 '-Dimethylbiphenyl 

Benzo[a]pyrene 
Figure 1 1 - Five CYP substrates of varying size. 
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ethylbenzene p-xyfene o-xyiene 




aa>dichlorotoluene 7-ethoxycoumarln 2-inethylnaphthalene 

Figure 13. Substrates used to parameterize aromatic oxidation. 
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-6J 

Predicted In Rate t-Butoxy 

Figure 14. Correlation between experimental reaction rates with t-butoxy radical and 

DFT energies. 
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